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I. SUMMARY 

In the following paragraphs, a brief summary of the advances in our understanding of 
the strengthening mechanisms of discontinuous composites will be given. 

All of the existing strengthening mechanisms that have been proposed to explain the 
increase in strength of discontinuous reinforcement composites are based on the assump¬ 
tion that the reinforcement does not deform. We have found that there is an increase in 
dislocation density within TiB 2 (the reinforcement in NiAl) as a result of deformation of 
the composite at 1273 K. I his increased dislocation density indicates that the TiB 2 is 
deforming during the deformation of the composite. 

It has also been observed that there is no apparent dislocation generation due to the 
difference in coefficient of thermal expansion between TiB 2 and NiAl. However, there is 
evidence of dislocation generation about the TiB 2 particles (due to a modulus difference! 
during plastic deformation. 

A finite element method investigation has demonstrated how and why the deforma¬ 
tion is localized in discontinuous reinforced composite. This local deformation has been 
observed in SiC/Al and TiB 2 /NiAl composite. Even though there is higher localized 
volume function of TiB 2 or SiC, this is the region (in the matrix) in which the higher 
density of slip lines and dislocations occur. 

The addition of discontinuous SiC to A1 alloys can result in a five-fold increase in the 
yield stress. The magnitud ;n the increase is obviously a function of the volume fraction 
and the particle size of the Sic. Previously, it was proposed that the strength increase due 
to SiC addition to A1 alloys was the result of change in the matrix strength, i.e. an increase 
in dislocation density and a reduction of subgrain size. The data obtained from a series of 
experiments indicate that dislocation density increases with an increase in volume fraction 
of SiC and decreases with an increase in particle size. The subgrain size decreases as the 
volume fraction increases and increases as the particle size increases. There is a good 
correlation between the microstructural changes in the matrix and the changes in the yield 
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stress of the composites. 
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Richard J. Arsenault 1 

Authorized Reprint 1988 from Journal of Composites Technology & Research, Winter 1 988 
Copyright American Society for Testing and Materials, 1916 Race Street, Philadelphia, PA 19103 

Relationship Between Strengthening Mechanisms and 
Fracture Toughness of Discontinuous SiC/AI Composites 


REFERENCE: Arsenault. R. J., “Relationship Between Strengthening 
Mechanisms and Fracture Toughness of Discontinuous SIC/AJ Com 

pushes, " Journal oj Composites Technology <4 Research. Vol. 10, No, 

4, Winter 1988, pp. 140-145. 

ABSTRACT: In all cates of discontinuous silicon carbide (SiC/AI) 
metal matrix it was shown that classical composite strengthening mech¬ 
anisms could not be used to explain the data. A high dislocation density 
model was proposed to account for the strengthening, but other factors, 
such as residual stress and texture, could contribute. These were inves¬ 
tigated. and it was shown that texture had no effect and residual stress 
reduced the tensile yield stress. The fracture process in the composites is 
controlled by fracture of the matrix up to SiC panicle sizes of 20 and 
for larger SiC size panicles fracture of SiC begins to dominate. The 
fracture process in the matrix is influenced by residual hydrostatic ten¬ 
sion and high density of dislocations generated at SiC/AI interfaces. 
Crack initiation fracture toughness does not depend on SiC particle 
size, for panicle sizes less than 20 pm. Crack growth resistance in¬ 
creases as the size of the SiC panicles increases. 

KEYWORDS: discontinuous SiC/ Al composites, strengthening mech¬ 
anisms, fracture toughness, residual stress, texture 

The framework of a mechanism is slowly evolving to account for 
the strengthening caused by the addition of silicon carbide (SiC) : 
to an aluminum (Al) alloy matrix. However, it should be kept in 
mind that the number of detailed investigations of these compos¬ 
ites is rather limited U-5). 

The basic strengthening mechanism is the high dislocation den¬ 
sity, which is produced as a result of the differences in the thermal 
coefficients of expansion between SiC and aluminum, and the 
small subgrain size that results. There could be other contributing 
causes to the strengthening such as 

(1) residual elastic stresses, 

(2) differences in texture, and 

(3) classical composite strengthening (load transfer). 

Therefore, a listing of the components of the strengthening 
mechanisms can be proposed 

A<Tyc = ACDISL "h i5sc i AOres A0TEX + A<7c0MP 

where Aaoisi is the increase in strengthening caused by the in¬ 
crease in dislocation density resulting from the differences in ther- 


'Profewor, Engineering Materials Group. University of Maryland. Col¬ 
lege Park, MD 20742. 

'In all cases the discussion is of discontinuous metal matrix composites. 


mal coefficients of expansion between the SiC and aluminum, Aa$c 
is the increase in strengthening caused by the reduced subgrain 
size. Acres ls the change in yield stress caused by presence of the 
thermal residual stress that results from the difference in thermal 
coefficients of expansion of the matrix and reinforcement. Aotex is 
the difference in strengthening caused by differences in text -e be¬ 
tween the 0 Vol% and the higher volume percent composites, and 
A°comp is the strengthening caused by classical composite 
strengthening, for example, load transfer. 

Dislocation generation about the SiC particle caused by differ¬ 
ences in the coefficient of thermal expansion was modeled by a 
simple prismatic punching model (6]. Figure 1 is a schematic of a 
particle that can take on the shape of a cube, plate, or whisker, 
from which several dislocations have been “punched" out. If the 
strengthening caused by this simple prismatic punching [6] is con¬ 
sidered then 20 Vol% spherical SiC with a diameter of 0.5 jim 
would produce a A a of 133.9 MPa, and 20 Vol% of whisker of l/d 
of 2 with a diameter of 0.5 /im (which is the approximate average 
diameter of the whiskers in the composite) would produce a A o of 
124.2 MPa. The model predicts a higher strengthening for the 
spherical morphology as compared to the whisker morphology- 
caused by the difference in volume of the SiC, not the volume frac¬ 
tion: the volume fraction of SiC remains the same. 

The change in strength caused by a large decrease in subgrain 
size can be obtained fTom the experimental investigations of 
McQueen and Jonas [7], 

An analytical model by Arsenault and Taya [4] based on an el¬ 
lipsoidal-shaped SiC particle in the aluminum matrix was devel¬ 
oped, which predicts that a tensile thermal residual stress should 
exist in the matrix for a whisker of l/d of 1.8 and that the longitudi¬ 
nal residual stress should be higher than the transverse residual 
stress. Also, in the case of composite containing spherical rein¬ 
forcement there should be zero net residual stress. 

A mechanism or mechanisms of fracture in these composites 
have not emerged as yet, and low ductility and fracture toughness 
of discontinuous SiC/AI composites remains a major obstacle to 
the practical application of these materials. Despite the significant 
improvement in the processing of SiC/AI composites in the recent 
years, fracture toughness is still in the range of 12 to 33 MPa • m 1 J 
\8-tl\. It is not clear how the interfacial SiC/AI bond strength, 
difference in CTE between SiC and aluminum, and size of the SiC 
reinforcement affect the toughness. It has already been established 
that the increase in the volume fraction of SiC panicles or whiskers 
or both and the increase of the matrix strength adversely influence 
the toughness of SiC/AI composites [10.12]- 

The purpose of this investigation was to consider the above men- 

© 1988 by the American Society tor Testing end Materials 
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FIG. 1— A schematic diagram of the particle and several prismatic 
punched dislocations. 


tioned three remaining strengthening items, that is, residual 
stress, texture, and composite strengthening. This paper will also 
examine the dependence of the toughness and tensile properties of 
the SiC/AI composites on size and morphology of the SiC 
reinforcement. 

Materials, Teat Methods, and Result* 

Two different SiC morphologies, whisker (//</ * 1.8 to 2.8 
where d =» 0.5 |im| and spherical (-0.5 nm in diameter), and 
three different matrices, 1100, 6061, and 7091 Al alloys, were 
used. In fracture toughness investigations, various sizes of SiC par¬ 
ticles were used, and they will be described later. 

The specimens used for the residual stress measurements were 
cylindrical. -12.5 mm in diameter by - 25 mm long, and had four 
flats 3 mm wide machined along the length of the cylinder. The 
specimens were then annealed for 12 h at 803 K and furnace 
cooled. The residual stresses were determined by measuring the 
changes in the lattice parameters for the various composites. There 
is a recognizable error in this method, but it was found that stan¬ 
dard techniques of using X-rays were completely useless. The stan¬ 
dard techniques assumed the specimen has a uniform tensile or 
compressive stress 1 13], along the length of the specimen. How¬ 
ever, in discontinuous composites this is not the case as Arsenault 
and Taya have shown [¥], the residual stress varies from point to 
point within the specimen. 

The actual experimental data obtained from an X-ray analysis is 
shown in Table 1. The X-ray results do indicate that the matrix is 
in tension and that the longitudinal residual stress is higher than 
the transverse residual stress. However, the X-ray data indicate a 


higher valu. of residual stress than predicted |4|. Also, the residual 
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TABLE l — The thermal residual stresses Itensilel as measured by X-ray 
diffraction techniques 


Material 

Transverse, MPa 

Longitudinal. MPa 

0 Vol% whisker SiC 

6061 matrix 

0.0 

0.0 

S Vol% whisker SiC 

6061 matrix 

34.5 

407.1 

20 Vol% whisker SiC 
6061 matrix 

S5.2 

227.70 

Wrought 

0.0 

0.0 


stress will effect the yield stress. The model of Arsenault and Taya 
(4) predicts that the yield stress in tension will be less than that in 
compression for whisker SiC reinforced composites. If the SiC is in 
the form of spheres there should be no difference in tensile or com¬ 
pressive yield stress. In the case of the spherical SiC reinforced 
composite there are tensile and compressive residual stresses in the 
matrix, but they cancel each other so the NET residual stress is 
zero. Table 2 lists experimental data obtained from 1100 alumi¬ 
num matrix specimens tested in the annealed condition. The sur¬ 
prising result is that in the case of spherical SiC composites, the 
tensile stress is higher than that of the compressive yield stress, 
whereas the difference predicted by the model is zero. However, 
the point to be made is that for the whisker case of > of. The 
model of Arsenault and Taya [4] successfully predicts the differ¬ 
ences in the tensile and compressive yield stress because of the 
thermal residual stresses. However, the model is completely incap¬ 
able of predicting the absolute magnitude of the increase in yield 
stress, for it does not have the capability of predicting the increase 
in the matrix strengthening. 

A "texture” [14] investigation was undertaken, to determine if 
the presence of the reinforcement would produce a difference in 
texture from that of the unreinforced matrix material. If there is a 
difference in the texture between the composite and unreinforced 
matrix material there will be a difference in the yield strength of 
the composite and the unreinforced matrix material [/5]. From a 
comparison of the data (Fig. 2) it is apparent that there is little 
difference in the texture of the 99.99% aluminum. 0 Vol%, 6061 
aluminum alloy and 20 Vol% 6061 aluminum alloy matrix com¬ 
posite. Therefore "texture" does not produce any difference in the 
yield strength of the composite and the unreinforced material. 

The remaining term to be considered is the Ao C omp- which is de¬ 
fined as classical composite strengthening or load transfer 
strengthening. The load transfer model developed by Taya and Ar¬ 
senault [5] as stated before predicts very little strengthening. This 
model, which is based on Eshbely’s model (76], handles load trans¬ 
fers in much greater detail than the shear lag model, and this cer¬ 
tainly takes into account the end effects [77]. For the case of whis- 


TABLE 1—Difference in yield stress. 


Matrix 

Vol% 

Ao,‘ MPa 

1100 

0 

3.8 

1100 

20. 

42.9 

1100 

20, 

-15.46 


‘Aa = af — a J 
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TEXTURE FOR EXTRUDED 8 ANNEALED 
20V% SIC 1100 MATRIX 

FIG. 2— The texture of composites containing various volume fraction 
of SiC 


ker composite with a I'd of 1.8 (which was obtained from 
h'stograms |/)) the model predicts that the composit" will have a 
slightly lower tensile yield stress than the matrix alloy with 0 Vol% 
SiC in it. The reason for this lower yield in the composite is due to 
the thermal residual stress. Figure 3 is a plot of yield stress versus 
volume fraction of SiC, and it is obvious that the experimentally 
observed strengthening is much greater than that predicted by the 
composite theory. Also, the composite theory predicts that if the 
shape of the reinforcement remains the same and volume fraction 
remains the same, the yield strength should be the same. However, 
experimentally, it wias found that by changing the particle size 
from - 5 to 250 |im, the yield stress decreases from - 1-72.5 to - 69 
MPa. 

Another investigation was undertaken to determine the effect of 
different matrices on the strength of the composite. Composite 
rods were obtained from ARCO Silag containing 0 and 20 Vol% 
SiC whisker with 3 different Al alloy matrices, 1100, 6061, and 
7091. The rods were machined into tension specimens and then 
annealed and furnace cooled. The difference in yield strength be¬ 
tween the 20 and 0 Vol% are listed in Table 3. There is a slight 
decrease in the magnitude of the increase in strengthening when 



FIG 3 —The yield stress as a function oj volume fraction of SiC. The 
solid line is of the experimental results The dashed line represents the yield 
stress as predicted by the shear lag model. 


TABLE 3 —Differences in yield stress for different 
matrices 


Matrix Material 



Annealed 

Aa." MPa 

% Atf, 

1100 

144 9 

420 

bOtil 

100 

181 

7091 

13! ! 

190 


"Ao, — a, jo v,,**, . 0, o vm-.. 


comparing 1100 and 6061 or 7091. It was expected that the higher 
yield stress of the 6061 and 7091 matrices, as compared to the 
1100. would result in a lower production of dislocations and a 
higher value of the thermal residual elastic stress. Both of these 
effects would reduce the increase in tensile yield stress of 6061 and 
7091 matrices, which is in agreement with experimental results. 

In the fracture investigation, specimens containing various size 
SiC particles were used. The SiC particles have a general shape of a 
platelets with a length to thickness ratio of 2 to 3. SiC particles of 
2.4, 3.2, 8, and 20 gm average length were mixed separately with 
1100 aluminum powder and hot pressed, extruded and hot rolled 
to produce tension test specimens and compact tension specimens 
(CTS). The composites were produced within the Metallurgical 
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Materials Laborato'y at the University of Maryland (/<$]. Also a 
SiC/AI composite containing 250-Mm-size SiC particles was pur¬ 
chased from DWA (DWA Composites Specialties. Chatsworth. 
CA). The volume fraction of the SiC reinforcement was constant in 
all composites and equal to 20 Vol%. An 1100 aluminum alloy was 
selected in order to minimize the influence of the alloying ele¬ 
ments. which would otherwise introduce additional complicating 
factors Tension samples were machined in the rolling direction 
and tested at a standard crosshead speed of 0.05 cm/min. The 
CTS were machined, with a T-L |/9) orientation, and tested using 
single specimen / integral testing technique in accordance with 
ASTM Test Method for a Measure of Fracture Toughness 
(E 813) 120). All specimens were tested in the annealed condition. 
It was initially assumed, based on previous results and theories [ 21] 
that the fracture toughness would increase with an increase in par¬ 
ticle size, thereby making it impossible to do valid K IC 1/9) mea¬ 
surements, since the sample thickness was confined to be 7 mm. 
Also, the energy separation technique developed at the University 
of Maryland (22) was utilized as a supplemental method of analyz¬ 
ing ioad-unload records to increase the confidence in the data ob¬ 
tained. The latter technique implies that the area under the load 
versus displacement curve, which corresponds to the work done by 
external load, can be separated into the stored elastic strain (po¬ 
tential) energy U„ the elastic energy U, released during crack ex¬ 
tension, and the plastic energy V p dissipated during crack exten¬ 
sion as shown in Fig. 4. 

The rate of the plastic energy dissipation / = 1 /B s dU„/da and 
the elastic energy release rate G - VB* dU,/da represent the 
plastic and elastic parts of the J integral, that is, J ~ I + C, where 
B n is the thickness of the CTS between the side grooves, and a is 
the crack length. Crack initiation fracture toughness can be deter¬ 
mined as K\c — (Cr£ c /I — v*) 1 2 (2J|. where E r is the composite 
Young's modulus and » is Poisson's ration (v = 0.J1). Also, a di¬ 
rect evaluation of Aiq is possible by taking maximum load P„ 
from the Ioad-unload record and substituting into the expression 
Kiq = P n / B N (W)' 2 /(a/*v) ( 14\, where w is the width of the CTS 
and values of/fo/w) are readily available. Crack growth fracture 
toughness is evaluated by dimensionless tearing modulus T. which 
is T = E c /al dJ/da, where o > is the composite yield stress and 
dJ/da is the slope of the stable crack extension portion of the J 
versus a plot constructed in accordance with ASTM E 813 (23.24). 
Crack extension was determined using the unloading compliance 
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FIG. 4— Schematic of the energy separation principle between two suc¬ 
cession unloadings. 
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technique f/2. /4). In order to verify the calculated crack extension 
values, the tested CTS were exposed to elevated temperature and 
then fractured. Physical crack lengths are within 5 to 10% from 
the values calculated by the compliance technique. 

Crack initiation fracture toughness measured as A, c and K, 0 is 
plotted as a function of the average SiC particle size and is shown 
on Fig. 5. 

The fact that both of these values show the same trend, that is. 
no dependence on the size of the SiC panicles, if the panicles are 
less than 20 jim in length, increases the confidence in the results 
obtained and supports the energy separation method as a new and 
powerful tool. The numerical difference (A| C = 18 MPa m 1 2 and 
K io s 23 MPa ■ m 1 2 ) between the values of K, c and K , 0 can be 
explained as follows. A finite notch with a root radius of « ISO *im 
was machined in CTS in order to stan the crack (as opposed to 
ASTM Test Method for Plane-Strain Fracture Toughness of Me¬ 
tallic Materials |E 399) requirement of a pre-fatigued crack starter 
resulting in a geometrically sharp crack [ 19]). Thus it is reasonable 
to expect higher values of A| 0 , which are obtained in accordance 
with ASTM E 399. Using the result of the study of the influence of 
notch acuity on the fracture initiation toughness |/0) in SiC/AI 
composites we may write 


/fio — Aici,i — K i 


(1 + r/2c? 2 
(1 + r/c) 


where Kj Cl ,> is the initiation toughness of the notch with radius r, 
K, c is the toughness of geometrically sharp crack, and c is the ad¬ 
justable constant related to the microstructure. Substituting values 
for Kior) = 23 MPa • m 1 2 , /(j c = 18 MPa • m 1 2 and r 3 150 pm, 
we obtain value for c S 20 pm. It is rather comparable with the 
size of the fracture process zone observed in 2A-pm average SiC 
particle size composite (25). Initiation fracture toughness for 250- 
pm SiC/Al composite is almost by the factor of 2 less than for the 
rest of the tested composites. This is apparently because of the pre¬ 
mature cracking of 250 -miti size SiC particles. 

Crack growth resistance measured as tearing modulus T and 
plastic part / of the J integral is plotted versus the average size of 
SiC particles in Fig. 6. 

The increase of the crack growth resistance with an increase of 
the SiC particle size means that there is an increase in energy dissi¬ 
pation during CTack extension in the composite with a larger size of 
SiC particles. As tension test data show the yield stress of the SiC/ 
Al composites drops with the increase of the SiC particle size. The 
same behavior was reported in literature for various alloys contain¬ 
ing second phase particles (26). The size of the plastic zone is in¬ 
versely proportional to the square of the yield stress (27). Thus as 
the size of SiC particles increase, the size of the plastic zone also 
increases resulting in the increase of dissipated plastic energy, 
which in turn increases the crack growth resistance. 


Discussion and Co Delusions 

If the strengthening components are considered again, it is now- 
necessary to consider the morphology of SiC. Table 4 is a listing of 
the strengthening components for the whisker case. The strength¬ 
ening caused by Ao-tex *nd Aocomp is equal to zero. The summa¬ 
tion of the predicted strengthening and the observed strengthening 
are in very good agreement. Table 5 is a listing of the strengthening 
components for the spherical SiC case. The strengthening caused 
by Aoqisl ntay be a little higher than the value given; however. 
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ARSENAULT ON SiC/Al COMPOSITES 


TABLE 4 — A luting of the strengthening components tor the whuker 
SiC 1100 aluminum matrix. 


A»tc = Aodisl + i Acres Aoiex * Aa t owp 

Aooisi = 11 7 
Aojc = 55 MPa 
Acres = 34.5 MPa 
Actex = 0 
Aocomr ~ 0 

Ao„ = 117 + 55 - 34.5 = 137.S MPa 

Aoyw e*r — 144.9 MPa 


TABLE 5 — A luting of the strengthening components for the spherical 
SiC in 1100 aiuminum matrix 


Aoyc — Aoqisl 4* Aosc i Acres 4 Aceex 4 Ao ( om, 
Acqise = 135 MPa 
Acsc = 55 MPa 
Acres = 0 
Acytx = 0 
Aocomp = 0 

Ao„ = 135 4 55 = 190 MPa 
Acys exr = 172.5 MPa 


again the agreement between the predicted and experimental 
results is very good. 

The data generated further support the concept that the 
strengthening mechanism is due to a higher dislocation density and 
a small subgrain size, which is the result of the difference in the 
thermal coefficient of expansion between SiC and aluminum. The 
thermal residual stress as measured by the X-ray technique indi¬ 
cates that the matrix is in tension. Also, a model based on load 
transfer is completely incapable of explaining the increase in 
strength due to the addition of SiC to the aluminum alloy matrix, 
in the composites used in the present investigation. The aspect ra¬ 
tio, in the whisker and platelet case, are too small to provide any 
effective local transfer. 

Since the fracture tougim.ss is independent of particle size and 
particle spacing (provided the particles are < 20 pm), the data sup¬ 
ports the contention that fracture process is controlled entirely 
within the matrix. In other words the fracture process (for the case 
of 20 pm and smaller SiC reinforcement composites) does not oc¬ 
cur by the nudeation and growth of voids of SiC particles. This 
does not mean that the reinforcement has no effect on the fracture 
process. The reinforcement can produce a high dislocation within 
the matrix; therefore the fracture toughness of composite is related 
or equivalent to the fracture toughness of the cold worked matrix 
alloy. Also, the reinforcement results in a hydrostatic tension in the 
matrix, which also further reduces the fracture toughness [25]. 
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EFFECT OF PARTICLE SIZE ON FRACTURE TOUGHNESS OF 
SIC/AL COMPOSITE MATERIAL"*" 

Y. Flom* and R. J. Arsenault 
Metallurgical Materials Laboratory 
Engineering Materials Group 
University of Maryland 
College Park, MD 20742 


ABSTRACT 

Discontinuous SiC/Al composites with SiC particles of different sizes were 

fabricated in order to study the role of particle size on the fracture process. 
The fracture process is confined to a very narrow band and takes place within 
the matrix in composites containing small SiC particle sizes. In the composite 
reinforced with SiC particles of 20 pn and above fracture of SiC begins to 
dominate. The matrix is influenced by the high density of dislocations 
generated at SiC/Al interfaces due to the difference in coefficient of thermal 
expansion (CTE) between SiC and the A1 matrix. Crack initiation fracture 
toughness does not denend on SiC particle size. Crack growth fracture toughrcss 
increases as the size of the SiC particle increase. 


I. INTRODUCTION 




During the past several decades there has been a considerable 
investigations of the fracture characteristics of two phase alloys 
second phase is a discreet precipitate or particulate^ The 


number of 
where the 
review by 


* Presently at Goddard Space Flight Center, NASA, Greenbelt, MD. 

+ This investigation was supported by the Office of Naval Research under Grant 
No. N00014-85-K-0007. 
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Shwalbe^*^ and more recent publications by Bates v 1J J , Gerberich 1 " 1 °' and Firrao 
and Roberti^^ are very informative on this subject. In ductile fracture, 
characterized by void nucleation and growth (VNG), the spacing between the void 
nucleating particles is generally considered to be a critical microstructural 
parameter that taken together with the tensile properties controls the toughness 
of a given material^^ This can be represented by the following simplified 
expression^ 


(15) 


,( 16 ) 


# 


_ r 


r 1/2 1/2 

!aa^e*Ef(N)J s 


( 1 ) 




♦ 


where is the plane strain fracture toughness, a is the numerical 

* 

coefficient, o is the yield stress, e f is the maximum strain acting at the 
crack tip, E is the Young's modulus, f(N) is some function of the strain- 
hardening exponent and s is the average inclusion spacing in the matrix. The 
experimental data collected by various investigators are usually given in terms 
of the volume fraction of the second phase particles (V) or less often in terms 
of s (see Fig. 1). 

The dimple morphology of the fracture surfaces is the most common 
observation made by a number of researchers studying the fracture process 
in discontinuous"*" SiC/Al composites^^ 32)^ This indicates that the VNG 
mechanism is active in these composites and therefore, corresponding theories 
(generalized by Eq. 1) should be applicable to describe the fracture process in 
SiC/Al composites. In the past, the investigations of the fracture toughness of 
SiC/Al composites were concerned with the matrix properties (heat treatment, 
amount of work, alloying elements, crack tip acuity) and the morphology of the 


+ In all cases we shall he referring to discontinuous SiC/Al composites. 
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SiC reinforcement (shape and volume fraction). 

The adverse influence of an increase of the volume fraction of SiC 
reinforcement on fracture toughness can be established rather well when one 
combines the data on fracture toughness testing available in the literature, as 
shown in Fig. 2. 

A rather strict geometric consideration requires that all three 

microstructural variables V, s and particle size d be related. There are 

Cl 1 1\ 33) 

several expressions relating V, s and d available in the literature'' ’ ’ 

Since there is no substantial difference between these expressions, we have 

(14) 

selected one derived by LeRoy et al v ' for the equiaxial particles: 

s = 0.77 dV _1/2 (2) 

The data in Fig. 2 can now be considered as a function of the average center-to- 

center SiC particle spacing s, providing that the average SiC particle size 

remains constant which is generally the case in commercial composites. There 

would be approximately a of factor 3 increase in if the data in Fig. 2 is 

employed along with Eq. 2, and assuming that the average SiC particle size is 

( 21 ) 

about 5 urn . Kotice that in this case we are changing s by varying V. 

Mathematically, a similar result can be achieved by keeping V constant and 
changing the particle size d. A similar trend is obtained, there is a linear 
increase of as d increases, assuming that voids are nucleated at SiC 

particles. Therefore, it appears that the fracture toughness of SiC/Al 
composites can be improved when the size of the SiC particles is increased, 
providing that at least one of the two following assumptions is met: 

1) voids are nucleated by SiC particles and/or 

2) the response of the matrix to changes of s remains the same regardless 





14 


of whether this change is caused by varying the volume fraction or the 
® size of the SiC particles. 

The purpose of this work, therefore, was to determine the influence of the 
size of SiC particles on fracture toughness of SiC/Al composites and to examine 
® the applicabilitv of the aforementioned assumptions. 

II. MATERIALS 


# 




Inert gas atomized 1100 A1 alloy (Valimet H-3) was mixed in succession with 
2.4, 3.2, 8 and 20 pm average size SiC particles obtained from the Norton Co. 
An 1100 A1 alloy was selected in order to minimize the influence of the alloying 
elements which would otherwise introduce additional complicating factors. The 
size distribution of the as received SiC powders were obtained by filtering onto 
0.4 pm size Nuclenore filters, and photographing the filters In the SEM. Figure 
3 shows histograms of SiC particles size distribution for all sizes used in 
present investigation. The SiC particles had an approximate platelet 
morphology. Mixtures of A1 and SiC powders were hot pressed, hot extruded and 

hot rolled to produce tensile test and compact tension specimens. Details of 

(34) 

the fabrication process are described elsewhere . The volume fraction of the 
SiC reinforcement was constant in all composites and equal to 20 volume percent. 
Following the rolling process the densities of the specimens were determined 
using a Buoyancy method per ASTM C 693-74. Density measurements showed that 
composites were within 99% of theoretical density. 

Also, a SiC/Al composite containing 250 pm size SiC particles purchased 
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from DUA* in the form of a 25 mm thick plate was sliced into plates ~ 6.5 mm x 
12.5 mm x 304 mm. The "as rolled" and sliced plates were machined into compact 
tension specimens (CTS). Design of the CTS was based on ASTM standard E 813 and 
also on the latest AST" developments in the area of elastic-plastic fracture 
toughness test methods ^ ^ as shown in Fig. 4(a). All specimens were 
annealed at 803 K for 12 hours and furnace cooled prior to the mechanical 
testing. 

TIT. TEST METHODS 

The CTS were tested using the single specimen J-integral test method per 

ASTM standard E 813. Also, the energy separation technique (EST) was utilized 

(37) 

as an additional tool in the load-unload records analysis '. The EST enables 
one to obtain values of plastic, I, and elastic, G, contributions to the total 
value of J-integral. It implies that the area under the load versus displace¬ 
ment curve which corresponds to the work done by external force can be separated 
into the stored elastic strain (potential) energy, U g , the elastic energy, Ug, 
released during crack extension, and plastic energy, U^, dissipated during the 

crack extension, as schematically shown in Fig. 5. A more detailed description 

(37 ) 

of the EST is outside the scope of this paper and can be found elsewhere 
The rate of plastic energy dissipation, I = 1/B n dU^/da, and the elastic energy 
release rate, G = 1 /dU^/da, represent the plastic and elastic parts of the J- 
integral, i.e. J * I + G, where B n is the thickness of the CTS between the side 
grooves and a is the crack length. Crack initiation fracture toughness can be 
determined as^^ : 


* DWA Composites Specialities, Inc. 
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K 


IC 


" (' 


G E 


c_^l/2 


1 - 


2 J 


(3) 


where E c is the composite Young's modulus and v is the Poisson ratio (v = 

0.31). Also a direct evaluation of K-j-q is possible by taking maximum load P*. 

(39) 

from the load-unload record and substituting into the expression' : 


K 


IQ 


P f(a/w) 
m 


B (w) 
n 


1/2 


(M 


where K is the width of the CTS as shown in Fig. 4(a) and values of f(a/w) are 
readilv available. Crack growth fracture toughness is evaluated by a 
dimensionless tearing modulus, T, which is equal 


fc dJ 

2 da 
o 


(5) 


where a is the composite vield stress and dJ/da is the slope of the stable 

y 

crack extension portion of J versus a plot constructed in accordance with ASTM 
E 813. Crack extension was determined by using the unloading compliance 

nc T O N 

technique^ ’ . In order to verify the calculated crack extension values, the 

tested CTS were exposed to elevated temperature and then fractured (heat tinting 
method). Experimentallv measured crack lengths were found to be within 5 to 10% 
from the values calculated by a compliance technique. 

The interaction between the fracture path and SiC particles was analyzed by 
counting the number of SiC particles along a random path (RP) and the fracture 
path (FP), and by comparing SiC particle densities and size distributions. 
Particle size is represented by the longest dimension of the SiC particle. 
Random path particle count (RPPC) was performed by counting the particles along 
the perpendicular lines forming a square grid which was placed on the 
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micrographs of the metallographically polished cross sections of SiC/Al 
composites. Fracture path particle count (FPPC) was done by incorporating the 
SiC particles touching both matching flanks of the crack on the cross sections 
of the tested but not separated CTS. In addition, the upper limit number of SiC 
particles fractured by the crack could be reasonably estimated since both 
matching halves of the same crack were present on the micrograph. The latter 
was used to determine the nercent of SiC particles fractured by the crack path 
using the following ratio: 

ft of particles fractured by crack path 
FPPC 

This technimie is described in greater details in Ref. 34. 

IV. EXPERIMENTAL RESULTS 

A. SiC Particles - Matrix Interaction 

The experimental data indicated that the ratio of SiC densities measured 
along RP and FP does not depend on the SiC particle size, i.e. the ratio 
RPPC/FPPC * constant = 0.75 for 2.4, 3.2, 8 and 20 pm average SiC size 
composites. This means that one finds more SiC particles along FP than along 
any other direction. Thus, the crack is attracted towards the SiC particles. 
In the case of the 250 pm SiC size composite there were too few SiC at the 
fracture surface to arrive at reproducible values, but ratio was small. 

The significance of the specific value 0.75 is not fully understood. 

The occurrence of the attraction or deflection is determined by the sign of the 
triaxial stress state. In the case of a tensile triaxial stress state, the 
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( 41 ) 

crack front is attracted toward the particles . The presence of the tensile 

(42) 

triaxial stress state has been demonstrated by Shi and Arsenault 

The upper limit percentage of SiC particles fractured by the crack path 

remains around 8\ for 2.4, 3.2, and 8 pm average SiC particle size composites 

and increases to * 25% for 20 pm SiC/Al composites. However, in the case of 
250 pm almost all SiC particles were fractured on the fracture surface. This 

result can he treated on the basis of the critical flaw size in the SiC used in 

the Griffith fracture mechanics. The probability of finding a critical flaw 

(43) 

size in a small SiC particle is less than in the laree one 


B. Fractographv 

Three characteristics were observed in the study of fracture surfaces: 1) 
fracture surfaces have a dimple morphology, 2) there are two dimple populations: 
the first, which is associated with SiC particles, and increases its size as the 
size of the SiC particles increases, and a second type consisting of very small 
dimples 1 pm) located in the space between the SiC related dimples, 3) cleav¬ 
age of SiC particles becomes more evident as their size increases (see Fig. 6). 
The shape of the dimples is rather equiaxial. 

Despite the correlation between the coarse dimples and SiC particles on the 
fracture surfaces of SiC/Al composites, not all the SiC particles were 
completely exposed. This corresponds to the earlier observations of the 

fracture surfaces of SiC/Al composites' » ’ , and indicates good bonding 
between SiC and A1 matrix. 
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C. Metallography 

The examination of the polished cross sections of the tested tensile and 
CTS did not reveal the presence of the voids below the fracture surface (see 
Fig. 7). This means, that the fracture process is confined to an extremely 
narrow band and there is no apparent damage zone adjacent to the fracture 
surface which generally exists in the classical void nucleation and growth type 
fracture. 

It is possible to stop a J-integral test and cut and polish the sample as 
schematically shown in Fig. 4b. It is then possible to determine how the crack 
propagates through the composite. (It should be clearly stated that crack was 
examined in the unloaded condition.) 

There were only a few observations of "pull out" of the SiC particles from 
the A1 matrix, which indicates a good SiC/Al interfacial bonding. Voids were 
not observed in the vicinity of and/or ahead of the crack tip. Instead, a 
series of short cracks could be seen in the matrix in front of the continuous 
crack as shown in Fig. 8. Apparently, this cracking takes place in the matrix 
ahead of the crack tip and it is believed that crack propagation occurs by 
connecting these discontinuous microcracks. 

It seems that the short microcracks are associated with the clusters of SiC 
particles (see Fig. 8). Even though the distribution of the SiC particles is 
predominately homogeneous, on the microscale there are islands of high and low 
density of SiC particles. The clusters of the inclusions are considered to be 
the sites where the damage level reaches extreme values causing the fracture of 
the entire system^^. The degree of plastic constraint within the clusters 
could be much higher than in the rest of the matrix as a result of: 1) high 
dislocation density due to the difference in CTE^^»45), and 2) increase of the 
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tensile triaxial stress state during the deformation due to the plastic 
constraint in the matrix between the particles. These factors make SiC clusters 
favorable for crack initiation sites. 


D. Crack Initiation Fracture Toughness 


Crack initiation fracture toughness measured as and is plotted as a 

function of the average SiC particle size and shown in Fig. 9. The fact that 
both Kj r and show the same trend, i.e. no_ dependence on the size of the SiC 

particles, increases the confidence in the results obtained and supports the 
energy separation method as a new and powerful tool. 

The numerical difference (Kj C « 18 MPa • and * 23 MPa • m^ 2 ) 

between the values of and K^q can be explained as follows. A finite notch 

with a root radius of * 150 pm was machined in CTS in order to start the crack 
(in contrast with ASTM E 399 requirement of a fatigue precracking to form a 
geometrically sharp crack). Thus it is reasonable to expect higher values of 

Kjq which are obtained in accordance with ASTM E 399. Using the result of the 

(23) 

study of the influence of notch acuity on the fracture initiation toughness 
in SiC/Al composites we can write 


K IQ ■ K IC(r) 


bet 1 + M : 

(1 + r/c) 


where is the intiation toughness of a geometrically sharp crack and c is the 

adjustable constant related to the microstructure. Substituting values for 
K IC(r) = 22 Kjc = 18 MPa • m^ 2 and r * 150 pm, we obtain c = 20 

pm. This value is within the range of the microcracking ahead of the crack tip 


in a 2.4 pm average SiC particle size composite. 
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Initiation fracture toughness of a 250 yra SiC/Al composite is almost a 
factor of 2 less than the rest of the tested composites. This is apparently due 
to the premature cracking of 250 pm SiC particles. 

E. Crack Growth Fracture Toughness 

Crack growth fracture toughness defined as tearing modulus, T, and plastic 
part (I) of J-integral versus the average size of SiC particles is plotted in 
Fig. 10. The increase of the crack growth toughness with the increase of the 
SiC particle size means that more energv is dissipated during crack extension in 
the composite with a larger size of SiC particles. 

V. DISCUSSION' 

In this section we will address crack initiation and propagation toughness 
as measured by and tearing modulus T, respectively. Throughout this section 

the attempt will he made to show that the size of the SiC particles in S1C/A1 
composite does not have a strong influence on the mechanisms affecting the 
portion of the touchiness of SiC/Al composite, but that the particle size has an 
influence on the tearing modulus (T), The influence of high dislocation 
densities on Kj^. and T will also be considered. Finally a brief description of 
the separation mechanism will be given. 

A. Crack Initation Toughness 

The experimental results of this investigation have indicated that crack 
initiation fracture toughness is independent of the SiC particle size 
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(within the range of SiC sizes tested). In order to understand this result it 
is instructive to consider local energy dissipation mechanisms associated with 
the crack extension in SiC/Al composite. We are using the term "local" to 
emphasize the area in the vicinitv of the crack tip. The plastic deformation in 
the bulk of the SiC/Al specimen is not included since it is separated out by 
using the energy separation technique (EST). 

The elastic enerev release rate G consists of the several terns that can 
account for enerev dissipation. These terms can be identified as: 


G = y + Y 

+ y + 

Y„_„ 



(7) 

s V 

PLzone 

CTS 



where y is the energv consumed in 
s 

creation of 

the 

new 

surfaces, 

Y is the 

V 

energv consumed in the formation 

of the void 

sheet 

’ ^PLzone iS 

the energv 

dissipated into plastic deformation 

around the 

crack 

tip 

anJ Yts 

represents 

crack tip shielding effects described by Ritchie 

et al. 

(46) 

. It is 

rather well 


/ Ofl \ 

established that the y tern' ' is verv small in metallic systems compared 

s ' 

with y„ + v . Evaluation of y„_ is quite complicated and requires the 

PLzone v CTb 

consideration of the extrinsic toughening mechanisms affecting the conditions at 
the crack tip in SiC/Al svsten. Following Ritchie's nomenclature^^ these 
mechanisms can be divided into four groups: 1) crack deflection and 

meandering, 2) zone shielding, 3) contact shielding, 4) local plastic 
deformation and void formation. 

In the next few paragraphs we will examine how the change in SiC particle 
size may influence these mechanisms. Before we proceed, however, one important 
comment should he made. The term "toughening" mechanisms used in this 
manuscript is the original term emploved by Ritchie^^). In order to avoid 
possible confusion we decided to retain this terminology even though the 




23 


presence of these mechanisms in SiC/Al composites does not result in the 
improvement of the toughness of the composite over that of the matrix material. 

1. Crack Deflection and Meandering 

Crack deflection by SiC particles was observed in the present work. 

( 47) 

Cotterell and Rice showed that the change in the stress intensity at the 

crack tip is a function of the kink angle 0. No observable differences in the 
crack deflection angle were noticed in the present work in SiC/Al composites 
containing small or large SiC particles. In addition to the kink angle it is 
possible that the kink length (O is also an important factor^ 14 ^. 

If the length, l, of the deflected portion of the crack is considered then 
the following analvsis (not very rigorous but still useful) can he performed. 
Using the idealized geometry of the crack deflection depicted in Fig. 11 the 
total length L of the kinked crack can be expressed as: 

L = 2MI , (8) 

where N is the total number of particles along the crack path. Again from the 
geometry in Fig. 11 one can see that 

* = + rid) • (9) 

where the term < 1 /2d + nd) represents the amplitude of the crack deflection, 
and n is the coefficient determining the extent of the crack deflection by the 
SiC particle. It seems that n depends on the elastic properties of the silicon 
carbide and aluminum and can be considered independent of the size of SiC 
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particle. It is possible to express N as: 

( 10 ) 

the volume fraction of the SiC 
size of the SiC particle. 


As one can see Eq. 11 does not contain parameter d, i.e. L is size independent. 
Therefore the analysis above show that the influence of the crack deflection on 
the toughness of SiC/Al system is independent of SiC particle size. 

2. Zone shielding 

Some of the mechanisms in this group such as transformation toughening and 
crack field void formation are not active in SiC/Al system. Others like crack 
wake plasticity, and the effects of dislocation density will be addressed later. 

The mechanism of microcrack toughening will be considered here again in 
terms of the SiC particle size. The microcracking (consisting of both types, 
cleavage of SiC particles and matrix microcracking) of the composite in front of 
the crack tip reduces the stress intensity at the crack tip if it occurs in the 
correct orientation^®^. It seems, therefore, that when SiC size increases, the 
increased cleavage of the SiC particles should result in further reduction of 
the stress intensity at the crack tip, i.e., the toughness (K^) should 
increase. On the other hand, the larger the SiC particles the more the total 


R V 


N = 


where B is the thickness of the specimen, V is 
particles (remains constant) and d is the 
Substituting Eq^. 8 and 9 into Eq. 10 we obtain: 


B V 


t = 0 


d + 








25 


crack path is consumed by the cleaved SIC particles. This means that the larger 
portion of the composite fails now in a brittle, SiC-like manner. Clearly, 
these two events are w’orking in opposite directions. It can be speculated, 


therefore, 

based 

on these qualitative 

considerations 

that as 

SiC 

particle 

size 

increases 

these 

two opposing events 

compensate each 

other 

and 

their combined 

effect is 

such 

that the influence 

of microcracking on 

K IC 

is SiC 

size 


insensitive. It is realized, however, that quantitative analysis should be 
perform to check the validity of this speculation. 

3. Contact Shielding 

This type of extrinsic toughening (due to wedging, bridging, and sliding) 
would be more relevant to cyclic loading, and therefore, will not be considered 
here. 


A. Local Plastic Deformation and Void Formation 


(49) 

Davidson has shown that y + y_, 

v PLzone 

factors in absorbing the energy around the crack, 
can be done following the approach suggested by 
the work of forming the void sheet: 


could be the most important 
The evaluation of y 

v 

Davidson^^ for calculation of 


V 

v 



W m 2 
_1_ o 

2 E 



( 12 ) 


°l e l 

where W is the work per unit area of void sheet formation, VI = - 

v ' ’ o m 

m(Cj ) 

Oj and are stress and strain picked from the stress-strain curve at some 


convenient point; m = n + 1, where n is the strain hardening coefficient, b is 
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the depth of the dimple and h is the height of the volume affected by the dimple 
(lq) 

formation' ’. 

Since we observed basically two populations of dimples, Eq. 12 has to be 

modified to account for this. V’e can write: 


W m 


"v- K 0 ( K ,(§c)%['4T-(ir) 


i o b~ , W m b 0 



03) 

where subscripts 1 and 2 represent small (size does not change) and large 
dimples (size increases as SiC particle size increases), respectively. 
Coefficients Kj and ^ represent the fraction of the total fracture surface 
occupied by small or large dimples. In order to estimate Kj and ^ we need to 
involve the results of the SiC particle count. Recall that RPPC/FPPC is equal 
to 0.75, which means that 1.333 times more SiC particles can be found on the 
fracture surface than on any random plane cut through the SiC/Al composite. The 
fracturing of the SiC particles also has to be taken into account. Thus: 

K 2 = 0.27(1 - X ) 04) 

where factor 0.27 comes from the product of (volume fraction) x 1.333 or 
0.20 x 1.333 = 0.27 and x I s the fraction of SiC particles fractured by the 
crack front. Kj = 1 - 0.27 = 0.73. Since the systematic determination of the 
dimple geometry was not part of this work the accurate calculation of Eq. 13 is 
not possible. It is possible, however, to make an order of magnitude type 
estimate using the approximate average number of dimple dimensions obtained as a 
result of tedious observation of fracture surfaces of the subject SiC/Al 
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composite. The shape of the dimples is considered to be an ellipsoid with the 
minor to major axis ratios of 1/3 and 1/5 for small and large dimples, accord¬ 
ingly. By letting d/3 < < d (where d is the largest dimension of the SiC 

(49) 

particle) and 0.3 pm < hj < 1 un for lower and upper bound we can estimate 

W : 
v 

-4 

W = 0.9 x 10 MPa • m for 2.4 pm size SiC, lower bound 
v 

= 1 x 10 Mpa • m for 20 pm size SiC, upper bound 

As we can see W v (or y ) remains practically independent of SiC particle size 

_3 

and represents only a small fraction of G = 3.2 x 10 Mpa • m obtained 

experimentally in the present work. Thus, the work to form the void sheet makes 
a rather insignificant contribution to the energy absorbed during crack 

extension in SiC/Al composite. 

It seems, therefore, that most of the energy is absorbed within the plastic 

zone around the crack tip (f^, ). In order to estimate this energv one has 

PLzone 

to know the distribution of the plastic strain in a direction perpendicular to 
the direction of crack propagation^^. Preliminary data indicates that in the 
region immediately adjacent to the crack the plastic strain is large, as 
evidenced by very large dislocation density^®^. Also it was observed that the 
dislocation density was the same in the immediate region of fracture for two 
different volume fractions of SiC reinforcements^^ (SiC particle size was the 
same). Based on these results we can assume that as the particle size 
increases, i.e. an interparticle spacing increases (V = const.), the density of 
dislocation in the immediate region of the crack remains the same. Therefore 

^PLzone 


is independent of particle size. 
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We have shown that the y , y , , and _which are the main 

v PLzone CTS 

contributing factors to G and directly related to K IC (Eq. 3), are independent 
of particle size. Contrary to this, the classial VNG theory predicts that 
should increase as the particle size increases. 

B. Propagation Fracture Toughness 

The increase in crack growth fracture toughness (measured as I and T) with 


an increase 

of the SiC 

particle size 

can be 

treated 

on the basis 

of 

the 

relationship 

between 

the 

dislocation 

density 

and the 

spacing 

of 

the 

SiC 

particles demonstrated 

by 

Arsenault and 

Shl< 45 >. 

From their work 

one 

can 

see 


that the increase in the tearing modulus with an increase in particle size is 
the result of the increased plastic deformation that can be accommodated in the 
matrix as the interparticle spacing increases. 

C. Separation Process 

Based on the data collected in the present investigation, the crack 
propagation in the SiC/Al composite can be described as follows. The crack 
extension takes place after a certain amount of the damage in the form of 
discontinuous microcracks is accumulated ahead of the continuous crack tip. The 
crack propagates when these microcracks are connected by microvoid coalescence 
which explains the dimpled appearance of the fracture surface. The crack front 
is attracted towards the SiC particles due to tensile triaxial stress state 
present around the SiC particles as a result of differential thermal contraction 
of the SiC and A1 matrix on the cooling down from the fabrication temperature. 








VI. CONCLUSIONS 


Frora the data generated in this investigation we arrived at five general 
conclusions. 

1. An increase of SiC particle size did not improve fracture toughness of 

SiC/Al composite. was independent of SiC particle size up to 

20 pm average size. 

2. The crack growth toughness, i.e. the tearing modulus and plastic work, 
increased as the SiC particle size increased. 

3. No evidence of void nucleation at SiC particles was found. 

4. SiC/Al composite is a more complicated system than it appears: i) the sites 
of void nucleation were not well defined, and ii) the two alternative ways 
of increasing interparticle spacing did not produce the same results: 

• If the particle size was constant, and the volume fraction was increased 

( 18—23 3?1 

and spacing decreased, then Kj£ decreased^ ’ 1 . 

• If the particle size and spacing were increased, and volume fraction was 
held constant, then remained constant. 

5. It can be speculated that K^ depended mainly on the volume fraction of SiC 
particles. 
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List of Figures 

Fig. 1. Fracture toughness as a function of the "process" zone size. When VNG 
mechanism is active the process zone size is equated with the average 
spacing, s, between void nucleating particles (Gerberich, Ref. 16). 
Various forms of Eq. 1 are employed by Gerberich to predict the 
fracture toughness values (solid and dashed lines) and compare them 
with the experimental data reported by Shwalbe and Spitzig. 

Fig. 2. Fracture toughness as a function of SiC volume fraction for SiC/Al 
composites. The shortage of the experimental data on the subject 
resulted in rather small number of data points presented in this graph 
(Refs. 18-23, 32). 

Fig. 3. Size distribution histograms of as-received 2.4 yn (a), 3.2 pm (b), 

8 im (c) and 20 yn (d) SiC particles. 

Fig. 4. Compact tension specimen (CTS) geometry (a) and the area of 
metallographic analysis (b). 

Fig. 3. Schematic of the energy separation principle between two successive 
unloadings. 

Fig. 6. Fracture surface of 3.2 ym (a) and 20 ym (b) average size SiC/Al 
composites. 

Fig. 7. Polished cross sections of 2.4, 3.2, 8 and 20 ym particle size tested 
SiC/Al composites (a, b, c, and d, respectively). 

Fig. 8. Polished cross section showing crack tip region in 20 ym SiC particle 
size tested but not separate SiC/Al composite compact tension specimen. 

Fig. 9. Crack initiation fracture toughness of SiC/Al composites measured as 
(a) and K ^ q (6). 

Fig. 10. Crack growth fracture toughness of SiC/Al composites measured as I (a) 
and T (b). 

Fig. 11. Idealized geometry of the crack deflection by SiC particles. 




MML 110 



s.ftm 


Fracture toughness as a function of the "process" zone size. 
When VNG mechanism is active the process zone size is equated 
with the average spacing s between void nucleating particles 
(Gerberich, Ref. 16). Various forms of Eq. 1 are employed by 
Gerberich to predict the fracture toughness values (solid and 
dashed lines) and compare them with the experimental data 
reported by Shwalhe and Spitzig. 
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Fracture toughness as a function of SiC volume fraction for 
SiC/Al composites. The shortage of the experimental data on the 
subject resulted in rather small number of data points presented 
in this graph (Refs. 18-23, 32). 
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Fig. 5. Schematic of the energy separation principle between two successive 


unloadings. 
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Fig. 11. Idealized geometry of the crack deflection by SiC particles. 
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THE EFFECT OF EXTERNAL MECHANICAL LOADING ON THE CHANGE OF 
THE MATRIX RESIDUAL STRESSES IN SiC/Al COMPOSITES 

N. Shi and R.J. Arsenault 
Metallurgical Materials Laboratory 
University of Maryland 
College Park, Maryland 20742 USA 

ABSTRACT 

The change of the average matrix residual stress as a function of two 
particular kinds of simple loading history was investigated, i.e. 
tension and compression. The major factor which controls the pattern 
of the change is presented. 

INTRODUCTION 

Previous experimental investigations of the deformation behavior of 
the short fiber reinforced SiC/Al composites have indicated that the 
constitutive behavior of these composites is distinctly different 
under applied uniaxial compressive and tensile loadings. It is 
generally found that the compressive yield strength of these composi¬ 
tes is larger than that of the tensile yield strength while the 
Young's modulus for compression is smaller than that in tension 1 ' 2 . 
It has been pointed out by various researchers 1-4 that the differences 
in response of the SiC/Al composite material to tensile or compressive 
load is closely related to the matrix thermal residual stresses. The 
previous analytical and experimental results indicate that in whisker 
reinforced composites the average thermal residual stress in the 
loading direction (z) in the matrix is in tension 2 * 4 . However, the 
change of residual stresses which may be a principal factor that 
affects the flow properties of the composites as a function of applied 
load has not been reported in the literature. It is therefore the 
primary goal of this investigation to analytically evaluate such a 
change in the magnitude of the matrix residual stresses. Throughout 
the analysis, it is understood that the term "residual stress" is 
referring to the phrase "the average matrix residual stress in loading 
direction Z." 

FINITE ELEMENT (FEM) MODELING PROCEDURES 

In our previous research 5 , it was shown that a two-dimensional FEM 
model could qualitatively reproduce the changes in characteristic 
constitutive and flow properties of the whisker reinforced SiC/Al 
composites when the system is subject to either a uniaxial tensile or 
compressive loading. Based on the assumption that the change of the 
flow properties are largely controlled by the change in residual 
stresses, it is also conceivable that the change of residual stress 
can be qualitatively predicted. 

In the present residual stress analysis, the ADINA finite element 
analysis code was used. The parameters used in this simulation 
correspond to annealed 6061 aluminum reinforced by Sic whiskers. The 
samples were initially thermally treated by setting a zero stress 
state at 773 K and then cooling down to 293 K and loaded either in 
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tension or compression to determine the materials' response under the 
influence of the thermal residual stresses. Such a response is 
characterized by macroscopically generating stress-strain curves which 
can be compared with experimental curves 5 and microscopically analyz¬ 
ing the change of the residual stresses after the external load is 

• removed. 

The composites are assumed to be two-dimensional infinite array of 
hexagonal packed parallelpipeds imbedded into the matrix. If the 
volume fraction of whisker is constant, it has been shown 5 that a 
small variation in interparticle spacing will not greatly affect the 

• material's property. Therefore, in the subsequent analysis, the 
interparticle spacing was not selected as a critical parameter. 


In the current analysis, 20V% Sic whisker composites with whisker 
aspect ratio equal to 4 were investigated. A multi-constraint bound¬ 
ary condition was employed so that only a quadrant of the material 
• which contains the Sic whisker, as shown in Fig. 1, was considered to 
represent the whole infinite array. 



RESIDUAL STRESS ANALYSIS 


Fig. 1 Diagram of the compo¬ 
site unit cell for FEM calcula¬ 
tion where the shaded area 
represents a quadrant of the 
Sic whisker and the remainder 
is Al. 


The average matrix residual stress is obtained in such a way that the 
stress values generated by FEM were averaged over the entire volume of 
• the matrix. It was found that the average matrix residual stress 
along the tensile axes will drop initially and then will increase or 
decrease depending on the nature of the external loading(Fig. 2). 



Fig. 2 A prediction of the 
change in the average thermal 
residual stress as function of 
eternal strain. C.F. is for 
initially tested in compression 
and T.F. is for tested in ten¬ 
sion first. 
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Fig. 3 Schematic diagram of 
the thermal residual stress 
arrangement where compressive 
residual stress is at the end 
of the fiber and the rest of 
the matrix is under tension. 


and Fig. 3 shows a schematic of the thermal residual stress arrange¬ 
ment where at the whisker end the matrix is under compression while 
the rest of the matrix is under tension. 

THEORETICAL EVALUATION OF THE CHANGE OF RESIDUAL STRESS 
UPON EXTERNAL LOADING 

In this section, an analytical investigation was undertaken in order 
to obtain a physical understanding as to why there is an asymmetrical 
behavior due to plastic straining, i.e. compressive loading is not a 
mirror image of the tensile loading 1 ' 3 . It was intended to depict the 
actual main physical process through the simplest model possible. 

It is the authors 7 belief that the effect of the volume mismatch 
between the Sic particle and the A1 matrix should be responsible for 
the asymmetric behaviors of the composites. Therefore, a simple 
discontinuous composite model based on such an idea was developed. It 
is assumed that a piece of Sic whisker is embedded into A1 matrix, as 
shown in Fig. 1, except there is no bonding between Sic and A1 along 
the whisker axes. Thus, there would be no shear load transfer and the 
boundary between tension and the compression zone is assumed to be the 
extension of SiC/Al interface along the whisker axes. Moreover, the 
shear stress at the zone boundary is assumed to be negligible. Also, 
a multi-constraint type of boundary condition was adopted as in the 
FEM analysis. 

By following the same physical cooling and the subsequent mechanical 
loading process as described in the previous FEM section, the thermal 
residual stress can be determined as: 

For the case of tension first: 
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For the case of compression first: 
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In both Eq. 1 and Eq. 2, o Wr and o mR represent the residual stress in 
loading direction in the whisker and matrix, respectively. E w and Ejj, 
are the Young's modulus for the whisker and matrix, E m ^, is the work 
hardening rate for the matrix, a w and a m are the coefficients of 
thermal expansion for the whisker and matrix, a is the yield stress 
for the matrix, and E c is the all-elastic Young's modulus of the 
composite without thermal residual stress as opposed to the apparent 
composite Young's modulus in which part of the matrix started plasti¬ 
cally deforming right at the beginning of the loading 


E c = 


(1-V T ) E m E w 
F c E w +F w E m 


+ V T E m 


where V-j. is the volume fraction of the tension zone, and 


F c = 


v c v w 

v w +v c 
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V c and V w are the volume fractions of compression zone and Sic whis¬ 
ker, respectively. Furthermore, o mR is the thermal residual stress 
in the tension zone. When AT is large enough to have induced plastic¬ 
ity which is true for the SiC/Al system 6 ' 7 
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and the average matrix residual stress (a m ) will be 




( 6 ) 











Fig. 4 A comparison 
between FEM and ana¬ 
lytical results of 
average residual 
stress as function 
of total strain. 


Fig. 4 shows the schematic curve for Eq. (6) if Eqs. 1 and 2 are 
considered. 


DISCUSSION 


From a consideration of Fig. 4, it is clear that the discontinuous 
composite model gives good qualitative agreement with the FEM data on 
the slope of the change of the residual stress. To be more specific, 
not only that the average residual stress drop is evident, but the 
slope change can also be correctly predicted. To understand the 
underlying mechanism for the change in residual stress, it is helpful 
to define a term as average matrix effective stress (o m ) based on a 
weighted averaging scheme where an expression for the loading rate 
[da z /d p) 2 may be used as a weighting function, where a z is the stress 
component in the loading direction corresponding to applied load p. 
Reevaluating the FEM results in the previous section by this scheme, 
it can be found that the ratio of the average matrix residual stress 
over the average effective residual stress (0m/°m e ff) is greater than 

[(°m/°m ef f)I = l*55j. This simply means that in order to maintain the 

volume compatibility, the compression zone at the whisker end will 
respond faster than the rest of the matrix, i.e. higher stress 
(strain) rate. This local stress (strain) rate inhomogeneity will 
later alter the residual stress arrangement after unloading. 


Furthermore, it is worthwhile to point out that in the case of com¬ 
pression first, there is a sign change in slope do m /d after the 

composite starts yielding, which is predicted by both FEM and the 
discontinuous composite model. This is an indication that the inter¬ 
action between tension and the compression zone becomes totally 
different. Once the macro-yielding is observed, a fact which has to 
be considered is that the matrix material is deforming plastically 
throughout regardless of initial sign of the residual stress. It is 
also interesting to point out that while the tendency of residual 
stress change is agreeable with both FEM and discontinuous composite 
model, the actual magnitude of the residual stresses predicted by the 
analytical model is much smaller than that predicted by the FEM. It 
is speculated that the absence of interfacial bonding in the model is 
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the key reason for the discrepancy. This discrepancy may be solved by 
assuming a set of effective materials' constant to the matrix material 
near the interface (tension zone). 


CONCLUSION 

In summary, (1) the predicted change of the average residual stress is 
asymmetrical, i.e. regardless of sign of the loading, the magnitude of 
the average residual stress drops upon initial elastic deformation, 
and it will increase or decrease depending on whether the loading 
history is compressive or tensile, and (2) the change of the residual 
stress predicted by FEM is generally comparable to that predicted by 
the theoretical model which is based on the volume mismatch between 
the reinforcement and the matrix. This indicates that such volume 
incompatibility is the controlling factor to the change of the resid¬ 
ual stress. 
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Abstract 

Crack initiation fracture 
toughness does not depend on 
Sic particle size. Crack 
growth fracture toughness 
increases as the size of the 
Sic particle increases. The 
fracture process is confined 
to a very narrow band and 
takes place within the 
matrix. In the case of 
composite reinforced with SiC 
particles of 20 ^m and above 
fracture of Sic begins to 
dominate. The matrix is 
influenced by a triaxial 
tension in the matrix and a 
high density of dislocations 
generated at Sic/Al inter¬ 


faces due to the difference 
in coefficient of thermal 
expansion (CTE) between Sic 
and the A1 matrix. 

1. INTRODUCTION 
If ductile fracture is 
characterized by void nucle- 
ation and growth (VNG), then 
the spacing between the void 
nucleating particles is 
generally considered to be a 
critical microstructura1 
parameter that taken together 
with the tensile properties 
controls the toughness of a 
given materials( 1-7 ). This 
can be represented by the 
following simplified expres¬ 
sion ( 7 ) 
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(1) 
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K IC = [ao y £ *Ef(N)] 

where K IC is the plane strain 
fracture toughness, a is the 
numerical coefficient, a y is 
the yield stress, ef is the 
maximum strain acting at the 
crack tip, E is the Young's 
modulus, f(N) is some func¬ 
tion of the strain-hardening 
exponent and s is the average 
inclusion spacing in the 
matrix. The experimental 
data collected by various 
investigators are usually 
given in terms of the volume 
fraction of the second phase 
particles (V) or less often 
in terms of s (see Fig. 1). 

The dimple morphology of the 
fracture surfaces is the most 
common observation made by a 
number of researchers study¬ 
ing the fracture process in 
discontinuous* SiC/Al compo- 
sites( 8-22 ). This indicates 
that the VNG mechanism is 
active in these composites 
and therefore, corresponding 
theories (generalized by Eq. 
1) should be applicable to 
describe the fracture process 
in SiC/Al composites. In the 
past, the investigations of 
the fracture toughness of 
SiC/Al composites were 


* In all cases, we shall be 
referring to discontinuous 
SiC/Al composites. 



2 3 4 5 6 7 8 


Fig. 1 Fracture toughness as 
a function of the "process" 
zone size. Various forms of 
Eq. 1 are employed to predict 
the fracture toughness values 
(solid and dashed lines) and 
compare them with the exper¬ 
imental data. 


MML 107 



SiC VOCUME FRACTION, % 


Fig. 2 Fracture toughness as 
a function of SiC volume 
fraction for SiC/Al composi¬ 
tes (Refs. 8-13,22). 
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concerned with the matrix 
properties (heat treatment, 
amount of work, alloying 
elements, crack tip acuity) 
and the morphology of the Sic 
reinforcement (shape and 
volume fraction). 

The adverse influence of an 
increase of the volume 
fraction of Sic reinforcement 
on fracture toughness can be 
established rather well when 
one combines the data on 
fracture toughness testing 
available in the literature, 
as shown in Fig. 2. 

A rather strict geometric 
consideration requires that 
all three microstructural 
variables , V, s and particle 
size, d, be related. There 
are several expressions 
relating V, s and d available 
in the literature ( 4 > 23 > 24 ) . 
Sic/Al composites. 

Since there is no substantial 
difference between these 
expressions, we select one 
derived by LeRoy et al.( 4 ) 
for the equiaxial particles: 

s = 0.77 dV*^ (2) 


20 

IS 

N 

~E 
• 
a 

a 10 

6 

X 

5 

0 

5 6 •> 8 9 

AVERAGE SIC PARTICLE SPACING, (im 

Fig. 3 Fracture toughness as 
a function of the average Sic 
particle spacing s. The 
experimental data in Fig is 
replotted using Eq. 2 and 
keeping Sic particle size d 
constant, equal to 5 urn. 
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Data in Fig. 2 can now be 
replotted as a function of 
the average center-to-center 
SiC particle spacing, s, 
providing that the average 
SiC particle size remains 


Fig. 4 Prediction of frac¬ 
ture toughness as a function 
of the average SiC particle 
size d. Volume fraction of 
SiC is constant and equaled 
to 20% (arbitrarily chosen). 
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constant which is generally 
the case in commercial 
composites. Figure 3 shows 
the dependence of K IC on s 
which was plotted using the 
data in Fig. 2 and Eq. 2 
assuming that the average Sic 
particle size is about 5 
Aim^ 11 ). Notice that in this 
case we are changing s by 
varying V. 

Mathematically, a similar 
result can be achieved by 
keeping V constant and 
changing the particle size d. 
In accordance with Eq. 2, by 
increasing d, we proportion¬ 
ally increase s. Thus, a 
plot similar to one shown in 
Fig. 3 should result if the 
size of Sic particles is 
increased and their volume 
fraction is kept constant 
(see Fig. 4). A similar 
trend can be obtained by 
plotting K IC as a function of 
s employing Eq. 1 (which is 
based on VNG theories of 
ductile fracture) assuming 
that fracture toughness of 
SiC/Al composites can be 
improved when the size of the 
Sic particles is increased, 
providing that at least one 
of the two following assump¬ 
tions is met: 

1) Voids are nucleated by 
SiC particles and/or 

2) the response of the 
matrix to changes of s 
remains the 


regardless of whether 
this change is caused 
by varying the volume 
fraction or the size of 
the SiC particles. 

The purpose of this work, 
therefore, was to determine 
the influence of the size of 
SiC particles on fracture 
toughness of SiC/Al composi¬ 
tes and to examine the 
applicability of the afore¬ 
mentioned assumptions. 

2. MATERIALS 

Inert gas atomized 1100 A1 
alloy (Valimet H-3) was mixed 
in succession with 2.4, 3.2, 
8 and 20 /im average size Sic 
particles obtained from the 
Norton Company. An 1100 A1 
alloy was selected in order 
to minimize the influence of 
the alloying elements which 
would other wise introduce 
additional complicating 
factors. As received Sic 
powders were filtered onto 
0.4 /m size Nuclepore filters 
for size distribution, and 
distributions were fairly 
broad (25) . Mixtures of A1 
and SiC powders were hot 
pressed, hot extruded and hot 
rolled to produce tensile 
test and compact tension 
specimens. Details of the 
fabrication process are 
described elsewhere( 25 ). The 
volume fraction of the SiC 
reinforcement was constant in 
all composites and equal to 


same 
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Also, a SiC/Al composite 
containing 250 irn size Sic 
particles was purchased from 
DWA* in the form of a 2 5 mm 
thick plate was sliced into 
plates „ 6.5 mm x 12.5 mm x 
304 mm. The "as rolled" and 
sliced plates were machined 
into tensile test and compact 
tension specimens (CTS). 
Design of the CTS was based 
on ASTM E813 standard and 
also on the latest ASTM 
developments in the area of 
elastic-plastic fracture 
toughness test methods( 26 • 27 ) 
and shown in Fig. 5(a). 

All specimens were annealed 
at 803 K for 12 hours and 
furnace cooled prior to the 
mechanical testing. 

3. TEST METHODS 
The CTS were tested using 
single specimen J-integral 
test method per ASTM E813 
standard. Also, the energy 
separation technique (EST) 
was utilized as a additional 
tool in the load-unload 
records analysis. The latter 
technique is quite new and is 
currently under intensive 
study by the fracture mecha¬ 
nics community ( 28 ) . The EST 
enables one to obtain values 
of plastic, I, and elastic. 



-SEE 0ET4IL A 



Fig. 5 Compact tension 
specimen (CTS) geometry (a) 
and the area of metallo- 
graphic analysis (b). 


* 
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G, contributions to the total 
value of J-integral. It 
implies that the area under 
the load vs. displacement 
curve which corresponds to 
the work done by external 
force can be separated into 
the stored elastic strain 
(potential) energy, U s , the 
elastic energy, U E , released 
during crack during crack 
extension, and plastic 
energy, U p dissipated during 
the crack extension. A more 
detailed description of the 
EST is outside the scope of 
this paper and can be found 
elsewhere (~ 8 ) . The rate of 
plastic energy dissipation I 
= 1/B n dUp/da and the elastic 
energy release rate G = 1/B n 
dU E /da represent the plastic 
and elastic parts of the 
J-integral, i.e., J = I+G, 
where B n is the thickness of 
the CTS between the side 
grooves and a is the crack 
length. Crack growth frac¬ 
ture toughness is evaluated 
by a dimensionless tearing 

modulus, T, which is equal to 
(29,30) 


E c dJ 



where a y is the composite 
yield stress and dJ/da is the 
slope of the stable crack 
extension portion of J vs. a 
plot constructed in accor¬ 
dance with ASTM E813. Crack 


extension was determined by 
using the unloading com¬ 
pliance technique( 26 > 29 ). In 
order to verify the calcu¬ 
lated crack extension values, 
the tested CTS were exposed 
to elevated temperature and 
then fractured (heat tinting 
method). Experimentally 
measured crack lengths were 
found to be within 5 to 10% 
from the values calculated by 
a compliance technique. 

A mini-tensile stage was 
obtained so that a precracked 
CTS could be load and crack 
propagation observed on the 
optical microscope. The 
sample was cut as shown in 
Fig. 5b and then metallo- 
graphically polished. The 
area in front of the crack 
was photomapped and then the 
CTS was loaded and photo- 
mapped again. 

4. EXPERIMENTAL RESULTS 
4.1 Crack Extension 
The CTS were section-polished 
and loaded such that the 
crack propagation could be 
observed under load. 

In the case of the composite 
containing small size Sic 
particles, i.e. <20 ^m, there 
are three general observa¬ 
tions: 1) there is very 

little, if not zero crack tip 
blunting or opening prior to 
propagation; 2) there is very 
little indication of void 







formation at the particle- 
matrix interface, and 3) 
there is very little indica¬ 
tion that particles are 
fractured as a result of the 
crack tip advance. 

4.2 Crack Initiation Frac¬ 
ture Toughness 

Crack initiation fracture 
toughness measured as K IC and 
Kjq is plotted as a function 
of the average Sic particle 


size and 

shown in 

Fig. 

6 . 

The fact 

that both 

K IC 

and 

Kjq show 

the same 

trend, 

i . e. , no 

dependence 

on 

the 


size of the Sic particles, 
increases the confidence in 
the results obtained and 
supports the energy separa¬ 
tion method as a new and 
powerful tool. 

4 • 3 Crack Growth Fracture 

Toughness 

Crack growth fracture 
toughness defined as tearing 
modulus (T) and plastic part 
(I) of J-integral is plotted 
vs. the average size of SiC 
particles in Fig. 7. The 
increase of the crack growth 
toughness with the increase 
of the SiC particle size 
means that more energy is 
dissipated during crack 
extension in the composite 
with a larger size of SiC 
particles. 

5. DISCUSSION 

5 . 1 Crack Initiation 

Toughness Kjg 

The experimental results of 


this investigation have 
indicated that crack initi¬ 
ation fracture toughness K IC 
is independent of the Sic 
particle size (within the 
range of SiC sizes tested) . 
In order to understand this 
result it is instructive to 
consider local energy dissi¬ 
pation mechanisms associated 
with the crack extension in 
SiC/Al composite. We are 
using the term "local" to 
emphasize the area in the 
vicinity of the crack tip. 
The plastic deformation in 
the bulk of the SiC/Al 
specimen is not included 
since it is separated out by 
using the energy separation 
technique (EST). 

The elastic energy release 
rate G consists of the 
several terms that can 
account for energy dissipa¬ 
tion. These terms can be 
identified as: 


G 


7 + 7 +7 
S V 


PLzone 


+7 


CTS 


(4) 


where 7 ^ is the energy 
consumed in creation of the 
new surfaces, j v is the 
energy consumed in the 
formation of the void sheet, 
is the energy dissi- 
pated into plastic deforma¬ 
tion around the crack tip and 
7 cts represents crack tip 
shielding effects described 
by Ritchie et al.( 31 ). It is 
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Fig. 6. Crack ini¬ 
tiation fracture 
toughness of Sic/Al 
composites measured 
as K IC (a) and k to 
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2 fracture toughness 

—i 

(measured as I and T) 

l 

-j with an increase of the 
Sic particle size can be 

I 

treated on the basis of 
the relationship between 
the dislocation density 
J and the spacing of the 
Sic particles. The 
crack growth fracture 
toughness is related 
primarily to the extent 
J of plastic deformation 

I 

I 

Fig. 7 Crack growth 
1 fracture toughness 

of SiC/Al composites 
measured as I (a) 
and T (b). 
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rather well established that 

the 7 S term( 29 ) is very small 

in metallic systems compared 

with 7 + 7 v Evaluation 

PLzone v 

of 7 cts ^i^e complicated 
and requires the consider¬ 
ation of the extrinsic 
toughening mechanisms affect¬ 
ing the conditions at the 
crack tip in the SiC/Al 
system. Following Ritchie's 
nomenclature( 31 ) these 
mechanisms can be divided 
into four groups: 1 ) crack 

deflection and meandering, 2 ) 
zone shielding, 3) contact 
shielding, and 4) local 
plastic deformation and void 
formation. 

5.2 Propagation Fracture 
Toughness 

The increase in crack growth 
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at a distance away from the 
crack. The amount of plastic 
deformation that can occur 
within the matrix depends 
upon the amount of plastic 
deformation in the matrix 
(i.c., due tc iCTE) prior to 
crack propagation. If we 
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equate dislocation density 
and the amount of plastic 
deformation then it is 
possible to determine the 
amount of plastic deformation 
permissible in the matrix as 
a function of particle size, 
due to external deformation, 

i.e., as a result of tensile 
or compact tension testing. 
In other words, if there is a 
low dislocation density in 
the sample prior to testing 
then the amount of deforma¬ 
tion tion that the matrix can 
accumulate will be large. 
Arsenault and Shi( 32 ) have 
shown this, by using a simple 
prismatic punching model. 

6 . CONCLUSIONS 
From the data generated in 
this investigation we can 
arrive at four general 
conclusions. 

1) Increase of SiC par¬ 
ticle size did not 
improve Kj C fracture 
toughness of SiC/Al 
composite. K IC is 
independent of SiC 
particle size up to 20 
/jtr average size. 

2) No evidence of void 
nucleation at SiC 
particles was found. 

3) SiC/Al composite is a 
more complicated system 
than it appears: a) the 
sites of void nucle¬ 
ation are not well 
defined, and b) the 
two alternative ways of 


increasing interpar¬ 
ticle spacing, S, do 
not produce the same 
results. Vt Si and 
K IC A ( d = constant). 
dt Si but K IC does not 
change (v = constant). 

4. It can be speculated 
that K IC depends mainly 
on the volume fraction 
of SiC particles. 
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MICROSTRUCTURE OF TiB 2 /NiAl* 

L. WANG and R.J. ARSENAULT 

Metallurgical Materials Laboratory, University of Maryland, College Park, Maryland 20742 

ABSTRACT 

The addition of 20 volume percent titanium diboride in particulate form (1-3 qm) to 
nickel aluminide (TiB 2 /NiAl) results in a twofold increase in the high temperature strength of 
NiAl. There are at least two theories that have been proposed to account for the high 
temperature strength of discontinuous reinforced metal matrix composites. However, they 
can not be adequately used as a basis to explain the observed strengthening. 

An investigation was undertaken of NiAl, 10 V% TiBJNiAl and 20 V% TiBJNiAl in 
the annealed condition and after deformation, allowed to cool slowly. There is a low disloca¬ 
tion density in the annealed samples and the dislocation density did increase slightly as a result 
of deformation However, deformation did produce some intriguing dislocation arrangements; 
for example, it was found that there was a high dislocation density within tire TiB 2 in the 
deformed higher volume fraction composites and the dislocation density within NiAl matrix 
was not uniform 

1. Introduction 

1.1 Background 

The interest in intermetallic compounds for high temperature applications extends 

back in time for a number of years [1] and there has been a continued interest [2-26]. There 

are numerous obvious advantages for the consideration of intermetallic compounds, i.e. 

they have good oxidation resistance at high temperatures. However, they are not widely 

used in structural applications, but they are, e.g. NiAl, used as coatings on super alloys for 

oxidation resistance. The reason that NiAl is not used in high temperature structural 

applications is that its strength is at least a factor of 10 weaker than conventional super 

This research was supported by the Office of Naval Research under contract No. 
N00014-K-85-0007. 








alloys. For example, at 1300 K and a strain rate of 10' 5 sec' 1 , the stress required for 
UDIMET 115 [27] is ~ 200 MPa, whereas for NiAl, it is < 20 MPa [26]. 

Another deterrent to the use of NiAl in structural application is its lack of ductility, 
especially in the large grain size condition [13]. The fracture toughness as measured by 
Charpy impact testing is about 0.08 ft.-lbs. (0.1J) [28]. It was assumed that the fracture 
toughness of SiC measured by Charpy impact testing was low, but in that case, the fracture 
toughness is 1 ft.-lb. [29]. 

The lack of high temperature strength and the low temperature ductility can probably 
be improved by the addition of TiB 2 , i.e. to produce a discontinuous composite of 
TiB 2 /NiAl. The major advantage of TiB 2 is that it is inert with respect to the NiAl matrix, 
and in high temperature applications, the matrix should not react with the reinforcement. If 
a reaction takes place, it is possible that reinforcement will be consumed during the 
reaction. 

1.2. Mechanism or mechanisms of strengthening 

The high temperature creep of composites as an indication of the high temperature 
strength of composites have all been analyzed in terms of a classical continuum composite 
strengthening mechanism [30]. Simply, this means that there is a reduction of the stress in 
the matrix due to the load transfer to the reinforcement. However, if these mechanisms are 
applied to a discontinuous spherical reinforcement, the load transfer is very small, i.e. 5%, 
which would mean that the predicted decrease in the creep rate of the spherical reinforced 
composite would be very small compared to that of the matrix material. 

If one attempted to perform a more rigorous analysis by the use of the Eshelby 
method, there is a restriction which is presently placed on the Eshelby method which 
requires that the matrix has to be treated as a neutonian fluid, i.e. the stress exponent (n) 
would have to be one (1). The observed stress exponents are much greater than one, in the 
range of 5 to 20 [31]. 

Another approach would be to consider the dislocation mechanisms that have been 
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proposed to predict creep behavior. However, very few detailed mechanisms have been 
• proposed which take into account the presence of large second phase particles [31,32]. 

From a consideration of a dislocation climb over second phase particles, it is possible 
to obtain the following equation [31] for the creep rate: 
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where D is the diffusion coefficient for the matrix, v is the Poisson’s ratio, X is the interpar¬ 
ticle spacing, E is the Young’s modulus of the matrix, p m is the density of mobile disloca¬ 
tions, o e is the effective climb stress, b is the Burger’s vector, r is the radius of the particle 
and kT are the Boltzmann’s constant and the temperature, respectively. If it is assumed that 
the second phase is approximately spherical, then there is a simple relationship [33] 
between x and r. 
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where Vf is the volume fraction of the reinforcement. As one would intuitively guess, for a 
constant Vf, the \ /D is greater for a larger size reinforcement, and e /D decreases with an 
increase in the volume fraction of reinforcement. If the actual e vs. o e is determined for a 
composite of 20 V% and 0.5 ^m radius particles in the NiAl matrix, then Eq. 1 predicts a 
stress for a € of lO^sec' 1 of <6.9xl0' 3 MPa. This stress is 3 orders of magnitude 
smaller than the strength of the matrix. There are other possible dislocation models which 
have been proposed by Blum and Reppich [32], but the result is that if the Orowan stress is 




low for the composite, the effect of the reinforcement on the creep rate is small. 

There have been a few TEM investigations of dislocations overcoming large panicles 
at high temperatures; however, details of the rate controlling mechanisms are not given 
[33]. 

It can be concluded that there is a large difference in the predicted high temperature 
strengths, i.e. as manifested by the creep rates of discontinuous composites as predicted by 
the continuum mechanics and dislocation models. If we consider the preliminary data 
obtained, there is approximately a factor of two increase of the strength of composite (20 
V%) vs. the matrix. 

The purpose of this investigation was to determine the microstructures of TiB 2 /NiAl 
samples which were deformed at high temperatures, and from these determinations obtain 
some insight into the possible strengthening mechanism. 

2. Materials 

The materials used in this investigation were purchased from the Martin Marietta 
Corporation, and consisted of 0, 10,20 and 30 V% TiB 2 in NiAl. The material was ground 
into compression samples with a height to diameter (h/d) ratio of l'/,-2 and were compres¬ 
sion tested at 1143 K and 1273 K at various strain levels. Samples also were annealed at the 
test temperature for periods of time from 2 hours to 12 hours. 

3. Testing methods 

The samples were tested in compression at various strain rates ranging from 10' 4 to 
10' 2 sec' 1 at 1143 and 1273 K. 

After compression testing and/or after annealing, samples were cut into 3 mm discs 
with a thickness of 0.5 mm by an electrical discharge machine. After lightly grinding to 
remove the debris on the surface, the discs were dimple ground until the thickness reached 
0.07 ~ 0.1 mm in the center. The final electro-polishing was performed by using 5% 
perchloric acid ethanol solution at 253 K or below until a small hole appeared at the center 
of the sample. The applied voltage was 50 ~ 70 V. 









Also, an investigation was carried out both optically and by scanning electron 
microscopy (SEM) to determine the particle size and distribution. This was also followed 
by an X-ray powder diffraction investigation to determine the content of TiB 2 . 

The transmission electron microscopy (TEM) was carried out at the high voltage 
electron microscopy (HVEM) facility at Argonne National Laboratory. 

4. Experimental results and discussion 

The yield stress (c>y) which is only slightly less than the ultimate stress, was deter¬ 
mined as a function of volume fraction of TiB 2 and strain rate. As expected, a y increases 
with increasing volume fraction (Fig. 1) and increasing strain rate. However, the increase in 
Oy with increasing volume fraction of TiB 2 is greater at faster strain rates. 

In Figs. 2-5 are SEM micrographs of 0, 10, 20 and 30 V% TiB 2 /NiAl samples. There 
are particulates in the 0 V% which are not TiB 2 . These have been tentatively identified as 
a-Al 2 0 3 from an examination of TEM selected area diffraction patterns. The other dark 
spots in Fig. 1 are holes which are created by an etching effect. It can be seen that the 
volume fraction of the second phase can be as high as 10 V%. The dark particles in Figs. 3, 
4 and 5 are the TiB 2 . In Figs. 6-8 are X-ray powder diffraction patterns obtained from 0,20 
and 30 V% TiB 2 powder samples. The purpose of the X-ray investigation was to obtain 
some insight as to what were in the SEM micrographs, i.e. an attempt to determine the 
nature of the other particles. 

In Figs. 9 and 10 are 0 V% samples deformed 30% at 1143 K (870° C). It is evident 
that the second phase particles are acting as barriers to dislocation motion. Figures 11 and 
12 are again of 0 V% TiB 2 /NiAl matrix samples showing the formation of a network which 
is prevalent in these alloys. Also, it is evident that the dislocation structure is not planar, i.e. 
no well defined slip plane, and the dislocations have convoluted shapes as shown in Fig. 12. 

Also, a phenomenon which is only observed in 0 V% NiAl matrix samples, especially 
the ones which were deformed at higher strain rates, are bands or regions of very high 
dislocation density which we define as shear bands. Figure 13 is an example of these shear 
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bands in a sample which was deformed at 1143 K. These shear bands appear to go straight 
across the grain boundary without any apparent shifting of the grain boundary or of the 
shear band. 

If we now consider the samples containing TiB 2 particulate, several things become 
evident. Figure 15 is a low magnification TEM micrograph of a 20 V% TiB 2 /NiAl which 
was deformed 30%. There is a very high dislocation density in the area with higher than 
average volume fraction of TiB 2 , whereas in regions of lower than average TiB 2 concen¬ 
tration, there is a very low dislocation density. Also, it is evident in these pictures that 
within the TiB 2 , there is a large number of dislocations. 

This difference in dislocation density which depends upon the local volume fraction 
of TiB 2 is evident in Figs. 16 and 17 where there are two different volume fractions of 
TiB 2 . Figure 16 contains a very high dislocation density in the matrix, whereas in Fig. 17, 
there is a very low dislocation density in the matrix. 

Returning to a consideration of plastic deformation which occurs in the TiB 2 , there 
are several unique observations as shown in Figs. 18-19. First, the TiB 2 particles are not 
uniformly deformed in all cases. Secondly, in the TiB 2 particles which contain a variable 
dislocation density, there is a variable dislocation density in the surrounding matrix. There 
is an inverse relationship if a portion of a particle has a high dislocation density, then the 
dislocation density in the matrix immediately surrounding this portion of the particle has a 
low dislocation density. Whereas, if a portion of the particle has a low dislocation density, 
then the matrix immediately surrounding this portion of the particle has a high dislocation 
density. 

As stated in the introduction, the purpose of this investigation was to determine the 
changes in microstructure of TiB 2 /NiAl composites deformed at high temperatures for the 
purpose of obtaining insight into possible strengthening mechanisms in this composite. The 
results obtained do not appear to support any of the presently proposed mechanisms. The 
theories of Taya et al. [30] which are based on the load transfer concept require that there 





should be zero or very little plastic deformation of the reinforcement. This requirement 
does not seem to be satisfied due to high dislocations within the TiB 2 particles. The 
theories of dislocation climb based on climb over second phase particles are inadequate, for 
they predict strengthening which is three orders of magnitude too small [31,32]. Finally, the 
theories of composite strengthening proposed by Arsenault et al. [34-42] based on an 
increased dislocation density and reduced subgrain size cannot be applied, for the 
increase in the dislocation density is too small and the change in subgrain size is not large 
enough. Therefore, an entirely new approach to causes of strengthening will have to be 
considered. 

5. Conclusions 

From a consideration of the data, there are several conclusions which we can arrive at. 
First, in regard to the NiAl matrix: 

(1) There are other particles present within the matrix and it could be as high as 10 
V%, 

(2) The Burgers vector is a {100} type with no defined slip plane, 

(3) Deformation at low strain rates results in no large change in dislocation densi¬ 
ties; at higher strain rates, shear bonds are observed in the 0 V% TiB 2 /NiAl. 

Second, in regard to the TiB 2 NiAl composites: 

(1) There is heterogeneous distribution of dislocations in the TiB 2 /NiAl 
deformed samples, 

(2) TiB 2 are deformed, 

(3) At present, a strengthening mechanism is not obvious. 
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Fig. 1 

Fig. 2 
Fig. 3 
Fig. 4 
Fig-5 
Fig. 6 

Fig. 7 

Fig. 8 

Fig. 9 

Fig. 10 

Fig. 11 

Fig. 12 

Fig. 13 
Fig. 14 

Fig. 15 


The yield stress at 1000° C (1273 K) of TiB,/NiAJ composites as a function of 
strain rate and volume fraction of TiB 2 /NiAl sample. 

Scanning electron microscope micrograph of 0 V% TiB 2 /NiAl sample. 

Scanning electron microscope micrograph of 10 V% TiB 2 /NiAl sample. 

Scanning electron microscope micrograph of 20 V% TiB 2 /NiAl sample. 

Scanning electron microscope micrograph of 30 V% TiB 2 /NiAl sample. 

X-ray powder diffraction pattern obtained from 0 V% TiB 2 /NiAl powder 
sample. 

X-ray powder diffraction pattern obtained from 20 V% TiB 2 /NiAl powder 
sample. 

X-ray powder diffraction pattern obtained from 30 V% TiB 2 /NiAl powder 
sample. 

Transmission electron microscope micrograph of 0 V% TiB 2 /NiAl sample, 
30% deformed, most of the dislocation having Burger’s vector of a[100]. 

Transmission electron microscope micrograph of 0 V% TiB 2 /NiAl sample, 
30% deformed, most of the dislocation having Burger’s vector of a[100]. 

Transmission electron microscope micrograph of 0 V% TiB /NiAl sample, 
30% deformed. Networks formed by reactions of [100] + [010] = [110] 
and[010] + [001] = [Oil], 

Transmission electron microscope micrograph of 0 V% TiB /NiAl sample, 
30% deformed. Networks formed by reactions of [100] + [010] = [110] 
and[010] + [001] = [011]. 

Schematic illustration for convoluted shaped dislocations. 

Transmission electron microscope micrograph of 0 V% TiB 2 /NiAl sample, 
30% deformed. Shear bond structure. 

Transmission electron microscope micrograph of 20 V% TiB 2 /NiAl sample, 
30% deformed. 


Fig. 16 Transmission electron microscope micrograph of 20 V% TiB 2 /NiAl sample, 
30% deformed. Highly concentrated TiB 2 area. 

Fig. 17 Transmission electron microscope micrograph of 20 V% TiB 2 /NiAl sample, 
30% deformed. Low concentrated TiB 2 . 

Fig. 18 Transmission electron microscope micrograph of 20 B% TiB 2 /NiAl sample, 
30% deformed. Not uniformly deformedTiB 2 particle. Deformed matrix 
always in contact with the dislocation-free end of the TiB 2 particle. 







Transmission electron microscope micrograph of the same area as in Fig. 17, 
except deformed end of TiB 2 particle with a practically dislocation-free area 
surrounding matrix. 
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Fig. 2 Scanning electron microscope micrograph of 0 V% 
TiB 2 /NjAl sample. 



Fig. 3 


Scanning electron microscope micrograph of 10 V% 
TiB 2 /NiAl sample. 











Fig. 4 Scanning electron microscope micrograph of 20 V% 
TiB 2 /NiAl sample. 



Fig. 5 


Scanning electron microscope micrograph of 30 V% 
TiB 2 /NiAl sample. 
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Fig. 8 X-ray powder diffraction pattern obtained from 30 V% 
TiB„/NiAl powder sample. 













Fig. 9 


Transmission electron microscope micrograph of 0 V% 

TiB /NiAl sample, 30% deformed, most of the dislocation 
having Burger's vector of a[100]. 



Fia 10 Transmission electron microscope micrograph of 0 V* 

y ‘ T iB /NiAl sample, 30% deformed, most of the dislocation 

having Burger's vector of a(l00]. 
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Fig. 11 Transmission electron microscope micrograph of 0 V% 
TiB 2 /NiAl sample, 30% deformed. Networks formed by 
reactions of [100]+[010]=[110] and[010]+[001]=[011]. 



Fiq 12 Transmission electron microscope micrograph of 0 V% 
TiB /NiAl sample, 30% deformed. Networks formed by 
reactions of [100]+[010]=[110] and [010]+[001]=[Oil]• 

























0% deformed. 
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Fig. 16 


Transmission electron microscope micrograph of 20 V% 
TiB 2 /NiAl sample, 30% deformed. Highly concentrated 
TiB 2 area. 
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Fig. 17 Transmission electron microscope micrograph of 20 V% 
TiB 2 /NiAl sample, 30% deformed. Low concentrated TiB 








Fig. 18 Transmission electron microscope micrograph of 20 B% 

TiB 2 /NiAl sample, 30% deformed. Not uniformly deformed 
TiB 2 particle. Deformed matrix always in contact with 
the dislocation-free end of the TiB 2 particle. 



Fig. 19 Transmission electron microscope micrograph of the same 
area as in Fig. 17, except deformed end of TiB 2 particle 
with a practically dislocation-free area surrounding 
matrix. 
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Dislocations in Titanium Diboride* 

L. Wang and R.J. Arsenault 

Metallurgical Materials Laboratory, University of Maryland, College Park, MD 20742 

ABSTRACT 

In this investigation, it was shown that titanium diboride would plastically deform at 
temperatures of 1300 K and at a stress level much lower than predicted by previous 
experimental data. The crystal structure of titanium diboride is a topologically hexagonal 
structure which allows for a number of possible slip systems. Experimentally, it was 
observed that three slip directions were present; [0001], <2110> and <2111 >. The 
probable slip planes are (0001), {1010} and {1011}. 

§ 1. INTRODUCTION 

Titanium diboride (TiB 2 ) is an interesting material. It has a high temperature 
strength which is greater than other ceramics such as silicon carbide (SiC) and aluminum 
oxide (A1 2 0 3 ) (Ramberg and Williams 1987). Also, TiB 2 is a reinforcement candidate in 
high temperature composite for it is inert with respect to potential matrix materials. One of 
the composite s> is which has received some attention is TiB 2 in nickel aluminide 
(NiAl) (Wang and Arsenault 1989). It was found that the addition of 30 volume percent 
particulate TiB 2 resulted in a five-fold increase in the yield stress at 1300 K. The reason 
why this increase occurs is not evident. 

A microstructural investigation was undertaken of these composites to determine the 
effect of the addition of TiB 2 on the dislocation arrangement and density within the NiAl 
matrix. It was assumed that the TiB 2 did not deform at 1300 K, for the yield stress at 1300 
K which is obtained by a lower limit extrapolation would be 2000 MPa. The stress 
concentration in the TiB 2 that could occur due to load transfer would not be large enough 

This research was supported by the Office of Naval Research under contract No. 
N00014-85-K-0007. 
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to cause yielding of TiB 2 at 1300 K. However, it should be pointed out that the TiB 2 used 
by Ramberg and Williams (1987) contained titanium carbide (TiC) precipitates. 

In the investigation of dislocation configuration and density (Wang and Arsenault 
1989) in the NiAl matrix of TiB 2 /NiAl composites deformed at 1273 K, it was observed that 
there was a high density of dislocations within TiB 2 particles. Therefore, a question arose 
as to what type of dislocations are within the TiB 2 particles. There appears to be a large 
number of possibilities. 

Nakaro, Imura and Takenchi (1973) and Nakaro, Matsubar and Imura (1974) in an 
investigation of the high temperature hardness of TiB 2 concluded that the primary slip 
system was {1010} < 1210 >. However, Ramberg and Williams proposed two other systems 
{10l'0}[0001] and (0001) < 1010>. 

The reason that so many possibilities exist is related to the crystal structure of TiB 2 
which is schematically shown in fig. 1. The crystal structure could be defined as a topologi¬ 
cally hexagonal structure (Parthe 1967), in which the boron atoms form a net connection 
only. For the case of packing of touching spheres (Laves 1956), the ideal c/a = 1.07, and 
the c/a = 1.063 for TiB 2 . From an examination of fig. 1, it is evident that there is no 
well defined close packed plane, but the closest packed plane is (0001) and the shortest 
repeat distance in this plane is [1120]. 

The purpose of the present investigation was to determine if an increase in dislocation 
density did indeed occur upon deformation of the TiB 2 /NiAl composite, and to character¬ 
ize the various types of dislocation Burgers vectors that were assumed to be present in the 
TiB 2 particles within the deformed TiB 2 /NiAl composite. 

§ 2. MATERIALS AND EXPERIMENTAL PROCEDURES 

The TiB 2 /NiAl composite used in this investigation was purchased from Martin 
Marietta Corporation. The composite contained 10 or 20 volume percent TiB 2 and the 
TiB 2 particles were of various shapes and size, but on average, the largest dimension was 
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3 ^m. 



Prior to testing, the compression samples with a length to diameter ratio of 2 were 
annealed at 1273 and 1673 K for periods of time of 2-12 hours. Variation in annealing 
times had no effect on the strength of the composite. Surprisingly, annealing at 1673 K 
resulted in a slight increase in strength. They were tested at 1273 K at various strain rates 
(details of the testing are given elsewhere (Wang and Arsenault 1989). The transmission 
electron microscopy (TEM) foils were obtained from the deformed composite samples and 
examined at both 100 kV and 1 MeV. Further details of the foil preparation are given 
elsewhere (Wang and Arsenault 1989). Since it is assumed that there will be a large 
number of possible Burgers vectors, it is necessary to use a number of operating diffraction 
vectors. However, if the relative intensity of the operating diffraction vector is so low, then 
the dislocations would be in very low contrast. Therefore, all micrographs in this paper 
were taken under operating conditions where the relative intensity is large enough to show 
the difference between the residue contrast of invisible dislocations and the image of visible 
dislocations. 


• §3. RESULTS AND DISCUSSION 

The initial portion of the investigation was to determine if the dislocation density 
within the TiB 2 increased as a result of deformation of the composite. In annealed 

• samples, there is low dislocation density both in the matrix and TiB 2 , as shown in Fig. 2(a), 
except in the case where some of the particles are in contact with each other. A comparison 
of TEM micrographs of annealed samples and deformed samples clearly shows that density 

• of dislocation within the TiB 2 is significantly higher in the deformed composite sample, as 
shown in figures 2(a) and 2(b). 

It was found that straight edges of a TiB 2 particle were almost always parallel to the 

• (0001) plane. The dislocations labeled A and B do not lie on the plane of the foil, but lie on 
a plane which is parallel to the bottom edge of the particle. In other words, dislocations A 



93 


& B lie on the (0001) plane which is perpendicular to the plane of the foil. Under 
conditions of 0001 and 2021 operating diffractions, dislocations labeled A in fig. 3(a) are 
almost completely invisible. Therefore, the A-labeled dislocations have Burgers vectors 
of a[1210]. Dislocations labeled B in fig. 3(a) and 3(c) which are only invisible under the 
operation diffraction 0001 has a Burgers vector of a[2110] or a[l 120]. 

If the electron beam is exactly parallel to the basal plane, then the projection of 
dislocation on the basal plane should be straight lines. However, there are exceptions if 
climb occurs as indicated by the arrow in fig. 3(a). Also, it has been found that stacking 
faults exist on the basal plane in TiB 2 (Merso, Lynch and Vahldiek 1968) Again, if the 
electron beam is parallel to the basal plane as in fig. 4(a), the dislocation on the basal plane 
should project as a straight line labeled A. If the sample is tilted as shown in fig. 4(b), then 
stacking fault ribbons are evident. The presence of stacking faults is not universal, only a 
few TiB 2 particles contained stacking faults. 

It is suspected that prism slip also occurs and that prismatic punching into the TiB 2 
particles from the interface occurs, i.e. the dislocations labeled C in figs. 5(a) and 5(b). If 
different operation diffractions are used (1010 and 2110), then prismatically punched 
dislocations become invisible (figs. 5(c) and 5(d), in areas that the Burgers vector of 
prismatically punched dislocations is c[0001j. 

There are dislocations with other Burgers vectors within this particle. The disloca¬ 
tions labeled E in figs. 5(a) and 5(c) are out of contrast under the diffraction vector 0001 
[fig. 4(d)]. The dislocations which have Burgers vectors of < 2110 > obviously do not lie on 
the basal plane. Therefore, it is possible that slip directions of the < 2110 > type can occur 
on other slip planes. 

Besides dislocations with Burgers vectors of <2110> and [0001], there are still other 
dislocations with different Burgers vectors, for example, the dislocations labeled D in fig. 5. 
It was difficult experimentally to find extinction conditions for these dislocations. Every 
operating diffraction condition was employed which had a required intensity in the 









(1210)*, (0110)*, (0111)* and (1213)* reciprocal planes. Not one of these operating 
diffraction conditions could make dislocations labeled D invisible. It is possible to obtain 
the Burgers vector of dislocations labeled D by considering the dislocation reaction which 
occurs at the arrow in fig. 5(b), (c) and (d). Dislocations labeled A which have a Burgers 
vector [2110], reacted with dislocation labeled D which results in the short segment of 
dislocation indicated by the arrow in these pictures. 

This short segment of dislocations has a Burgers vector [0001], for it followed the 
same extinction conditions as the dislocation labeled C which has a Burgers vector of 
[0001]. Therefore, the reaction can be written as 

b A + b D = b c 

The subscripts denote the dislocations labeled in fig. 5. Since b^ is [2110] and b c is [0001], 
then the dislocations labeled D must have a Burgers vector of the type < 2111 >. The reason 
why it is so difficult to find extinction conditions for these is that only 0110 diffractions 
satisfy extinction conditions and the intensity requirement mentioned above. 

The possible combinations of slip directions and slip planes are listed in Table 1. The 
three directions were determined in this investigation, and the planes listed are the three 
of highest atomic packing. 

§4. CONCLUSIONS 

From the data obtained in this investigation, it is possible to arrive at the following 
conclusions. 

• The TiB 2 particles within a deformed TiB 2 /NiAl composite can be plastically 
deformed as evidenced by the increased dislocation density within the TiB 2 particles. 

• Three slip directions were observed within the TiB 2 particles; they are <2110>, 
[0001]and <2111 >. 

• Dislocations with Burgers vectors of the type < 2110 > on the basal plane could split 
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into two particles with a stacking fault between the particles. 

Most of the dislocations in TiB 2 had Burgers vectors of [0001] or <2111>. 
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Table 1 

Slip Systems in Titanium Diboride 



Slip Directions 
[ 0001 ] 

< 2110 > 

< 2111 > 
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Fig. 1 Ball model of TiB 2 . Boron atoms form a net connection only. Titanium atoms sit in 
the open center of the net. (Westbrook 1967). 
































Fig. 3(c) 

Fig. 3(a-c) TEM micrographs of a TiB 2 particle in deformed 20 V% TiB 2 /NiAl composite 
with different operating diffractions, (a) dislocation with Burger’s vector of a < 2110 > type 
on the basal plane are visible, (b) all a<2110> type dislocations are invisible, (c) only 
dislocations with Burger’s vector of a[1210] are invisible. 








Fig. 4(b) 


Fig. 4(a-b) TEM micrographs of a TiB 2 particle in deformed 20 V% TiB 2 /NiAl compo¬ 
site. Stacking faults on the basal plane, (a) basal plane is parallel to the electron beam, (b) 
basal plane is not. 
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Fig. 5(c) 


Fig. 5(d) 






Fig. 5(a-d) TEM micrographs of a TiB 2 particle in deformed 20 V% TiB 2 /NiAl compo¬ 
site, (a) many beams condition. All the dislocations are visible, (b) all of the < 2110 > 
type dislocations are invisible, (c) and (d) prismatically punched dislocation out of contrast. 
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Abstract 

Tensile tests of SiC/6061 Al composites containing various volume fractions of whiskers 
or particulates (20, 5 and 0V%) showed that for high volume fraction samples (20 V%), the 
fracture process was very localized, i.e. a very narrow neck. As the volume fraction of 
whiskers or particles decreased, the deformed region spread out. One might expect that the 
microstructure should correspond to the macroscale changes. In the highly deformed region 
the dislocation density is expected to be higher, in the less deformed regions the dislocation 
density should be lower, and if the deformation is very localized, then the high dislocation 
density should also be limited to a very narrow region 

Overall, there is good agreement between the microstructure (dislocation density) change 
and the macroscale deformation of SiC/Al composite tensile samples. The mechanism 
proposal to account for this change in deformation behavior as a function of volume fraction 
of SiC in Al is related to the expansion of the plastic zone (due to differences in thermal 
coefficients of expansion between SiC and Al) when the external stress is applied. Also, the 
localized deformation is related to localized clusters of SiC particles. There is a cooperative 
effect which leads to a region of very localized plastic deformation. 

1. Introduction 

There are many unusual deformation characteristics of SiC/Al discontinuous metal 
matrix composites. Also, the ductility and fracture toughness of SiC/Al composites are very 
low [1]. Generally, examinations of the fracture surface do not reveal the presence of SiC 
particulates or whiskers in the quantities corresponding to their volume fraction [1]. The 
fracture process is very localized, very little "necking" is observed (in high volume fraction 
composites), but the fracture surface has both brittle and ductile characteristics. There is 

This research was supported by the Office of Naval Research under Contract No. 
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no indication of fracture of SiC (when the size of SiC is less than 10 /im). The fracture may 
initiate at large intermetallic inclusions, clusters of SiC particulates or whiskers and voids 
associated with the SiC particulate or whiskers [1]. 

The results [2] of the tensile tests of SiC/6061 A1 composites with various volume 
fractions of whiskers or particulates (20, 5 and 0 V%) indicated that for high volume 
fraction samples (20 V%) the fracture process was very localized. As the volume fraction of 
whiskers or particulates decreased, the deformed region spread out. Figure 1 demonstrates 
how "necking" changes with a decrease in volume fraction. In the case of the 20 V% 
composite the reduction in area at the fracture surface was only 3%, and 5 mm away from 
the fracture surface the reduction in area was zero. However, in the case of 5 V% composi- 



..... , tes at the fracture surface, 
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the reduction in area was 
14% and 5 mm away from 
UNDEFORMED the fracture surface, the 
reduction in area was 4%. 



Fig. 1 Reduction in area at various distances from 
fracture surface for SiC/Al composite tensile samples. 


One might expect in a 
highly deformed region, the 
dislocation density is 
expected to be higher and in 
the less deformed region 
which is very localized, the 
dislocation density should be 
lower. In other words, the 
high dislocation density 
should be limited in a very 
narrow region. It is also of 
interest to know that 











previous work [1] has indicated that along the fracture path, the local density of SiC 
particles is higher than the average volume fraction of the SiC. In other words, the fracture 
path is attracted to the particle clusters. However, it is still unclear whether such a "particle 
friendliness" is only restricted tc the fracture process or if it also governs the onset of plastic 
deformation. 

To investigate the underlying mechanism for such a phenomenon, a set of two- 
dimensional finite element method (FEM) meshes were generated to simulate the effect of 
inhomogeneous distribution of the SiC reinforcement whisker; that is, the effect of 
clustering. Furthermore, it is also of importance to know the plastic deformation process 
around each SiC whisker. Therefore, a three-dimensional FEM mesh was generated to 
monitor the behavior of plastic zones around whiskers. 

There have been numerous investigations and several theories of the ductile fracture 
of two phase alloys where the second phase is a discrete-non-deforming precipitate or 
particle [3-19], All of the theories are based on the assumption that a void nucleates at the 
interface of the non-deforming particle and then grows, leading to ultimate fracture. An 
examination of transfer sections by optical and scanning electron microcopy methods 
indicated there are very few voids near the fracture surface. It is possible that the voids 
nucleated are very small and can only be observed by transmission electron microscopy 
techniques. 

To determine whether this was a valid assumption, an investigation was proposed to 
determine the dislocation density along the tensile axis of the gauge section by using the 
transmission electron microscopy to connect the relation of macroscale and microstructure 
change to the plastic deformation. Also, the purpose of this investigation was to determine 
if there was a higher density of voids near the fracture surface in SiC/Al composites. 

2. Experimental Procedures 
2.1. Materials 

SiC whisker or particulate reinforced 6061 A1 alloys composites were used for this 









investigation. The materials were obtained commercially. 

The composites with 20, 5 and 0 V% SiC whiskers (SiC^/Al) were purchased from 
ARCO SILAG (presently, Advanced Composite Materials) in the form of extruded rod 
15.5 mm in diameter. The composites with 20 and 5 V% SiC particulates (SiCp/Al) were 
purchased from DWA Composite Specialties in the form of a plate 15.9 mm thick. All 
materials were fabricated by the powder metallurgy method. The whiskers [3] were small in 
diameter, ranging from 0.1 to 1 pm, and the average aspect ratio was 2. The particulates 
were platelet in shape, 5 to 7 pm long and had an aspect ratio of 2 to 3. 

Sample blanks were cut from the rods and plates, then machined into tensile samples. 
The doubly reduced samples were required to ensure that the fracture occurred in the 
gauge section and not in the transition region [20]. The dimensions of the samples are 



Fig. 2 Dimensions of a round tensile sample (in millimeters). 


given in Fig. 2. Those that were to be heat treated were given a modified T6 heat treat¬ 
ment, which means a solutionizing step of three hours (modified) as opposed to the 
standard one hour. The modified heat treatment was used in case a divorced eutectic 
formed during the processing. Other were annealed for 12 hours at a solutionizing 
temperature of 810 K and furnace cooled. 

2.2 Tensile testing 

All samples were tested in an Instron testing machine at room temperature. Most of 










the samples were tested to failure, however, some were tested to a strain where a measur¬ 
able "neck" occurred in the sample. 

2.3 TEM thin foil preparation and examination 

After the tensile test, a 3 mm diameter cylinder was electrical-discharge-machined 
and cut longitudinally (parallel to the tensile axis) out from the gauge section. Then several 
slices of a thickness of ~ 0.6 mm disk were cut transversely (along the tensile axis) from the 
cylinder by a low speed diamond saw. Each disk was fixed on a brass block by DUCO 
cement, then mechanically thinned on a rotating water flooded wheel covered with 600 grit 
paper to remove the diamond saw cutting damage and reduce the thickness to approxi¬ 
mately 0.12 mm. The reason for using DUCO cement rather than wax (the ordinary way) 
was to avoid heating the 3 mm disk during the mounting and dismounting of disks from the 
holders. Again by using DUCO cement, the disk was fixed on a holder and then the disk 
was dimpled by using diamond paste down to one micron in size forming a disk with a 
thickness of about 0.025 mm at the center of the dimpled region. The final thinning was 
carried out by using argon ion plasma bombardment operating at five kV and an ion 
current of 25 microamperes and a disk inclination of 11 degrees to the ion beam. 

All thin foils were observed in a AEI EM7 type high voltage electron microscope 
operated at 1 MeV at the Argonne National Laboratory. Some of the thin foils were also 
observed in a JEOL 100 CX electron microscope at 100 kV. The liquid nitrogen cold trap 
was used to minimize contamination. The double tilt stages were used in both microscopes. 
However, both microscopes are not equipped with an eucentric stage. 

There are several factors which will influence the accuracy of the value of the 
dislocation density, such as difficulties in determining the foil thickness. Other than this 
difficulty, one major problem in this part of the investigation was that, in case of very high 
dislocation density (>10 16 m' 2 ), the measurement will be less accurate because of the 
uncertainty of the real positions of dislocations on the TEM micrographs. 
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2.4 Dislocation density calculation 

The dislocation densities were determined by the method proposed by Keh [21]. In 
this method, a grid of two sets of lines normal to each other, but with different spacing, is 
placed over the TEM micrograph and then the intersections of dislocation if given by 

a = (Nj/Lj + N 2 /L 2 )t (1) 

where L x and L 2 are the length of the grid lines normal to each other, N x and N 2 are the 
intersections along the lines and t is the thickness of the area where the micrograph was 
taken. 

All the micrographs were taken from the area away from the foil edge to avoid 
underestimating the dislocation density. Generally, three to four different locations in each 
foil were examined, and micrographs were taken from three to five different areas in each 
location. 

The upper limit [22] of dislocation density calculated by using Eq. (1) is about 10 15 ~ 
10 16 m' 2 . In some micrographs, dislocations were so dense, it was impossible to see the 
individual dislocation. The dislocation density was estimated by using the width of the 
micrograph divided by width of the dislocation image. Clearly, the value of dislocation 
density was underestimated. 

2.5 Slip line analysis 

From melted 250 pm SiC particulate reinforced composite, rectangular samples 8 mm 
wide by 50 mm long and 1.5 mm thick were electrical-discharged-machined and polished. 
The samples were glued to a E. Fullam tensile substage model No. 18202 JEOL JSM-35. 
The substage was mounted onto a specially constructed holder to fit on to a Zeiss ICM-405. 
The entire effective gauge was photographed. Then the sample was deformed a small 
amount and the area which contained slip lines was re-photographed, and again the sample 
was deformed and re-photographed where slip lines occurred. This process was repeated 









no 

until the sample fractured. 

2.6 FEMprocedures 

2.6.1. T wo-dimensional FEM procedures 

In the present study, the composites are represented by a periodic array of trans¬ 
versely aligned whiskers as shown in Fig. 3, i.e. an infinite two-dimensional array of SiC 
whisker embedded in the aluminum matrix in plane strain condition. In order to consider 
the effect of clustering, the array is transformed in such a manner that an array of periodic 
geometric cluster centers was constructed and the neighboring whiskers (in this case, four 
whiskers [Fig. 4]) were attracted toward their cluster centers so that a periodic array of 
clusters were formed as shown in Fig. 5. The degree of clustering was defined by the value 
of Sy/V and S^/H in Fig. 5. 

Applying above geometric transformation scheme, a commercially available ADINA 
finite element materials (FEM) code was employed. To represent the reality, the material 
is assumed to be stress-free at annealing temperature, i.e. 773 K, and is subsequently cooled 
to room temperature. Therefore, the effect of thermal residual stresses is implicitly 
considered. After thermal treatment, the composites, clustered and unclustered, were 
uniaxially loaded in whisker directions. Due to the periodic arrangement, the condition of 
symmetry can be so utilized that a "unit cell" was selected as a basic building block such that 
only a quadrant of whisker (cluster) is contained as shown in Figs. 6(a) and (b) , and 
multi-constraint boundary conditions were employed to represent the symmetry conditions 
at the boundary. 

2.6.2. Three-dimensional FEM procedure 

In this case, the reinforcement particles are hexagonally distributed through three- 
dimensional space. Because of the symmetry, the response of the composite can be fully 
represented simply by considering a building block which only contains one-eighth of the 
reinforcement particle, as shown in Fig. 7 and, in consideration of symmetry, multipoint 
constraint equations were employed in the boundary conditions. As in the case of two- 
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Fig. 3 Periodic hexangonal arrangement of the SiC reinforcement of the whisker. The 
shaded areas are SiC whiskers. 



Fig. 4 This figure shows how SiC whiskers move to their geometric cluster center to form a 
periodic array of clusters. Whiskers with dotted boundaries represent the final locations of 
the whiskers after the movement indicated by the arrows in the figure, and Sy and Sjj show 
the degree by vertical and horizontal clustering. 















































































Fig. 6(b) A unit cell for the clustered whisker configuration employed by the FEM analysis 
where Sy and are defined the same way as in Fig. 4. 



Fig. 7 A unit cell employed for three-dimensional analysis. 










dimensional analysis, the composites were assumed stress-free at 773 K which is the 
annealing temperature for SiC/Al composites, and then the sample was subsequently cooled 
down to room temperature. Due to a large difference in thermal expansion coefficient 
between SiC and Al, severe thermal residual stress would be generated around particles. 
As a result, it is apparent that a plastic zone will form around each particle. Upon further 
uniaxial loading, such a plastic zone will certainly respond to the external load. Depending 
on how such a response would be, it will be helpful to understand the deformation 
mechanism of the composites. Therefore, it is our intention to monitor the geometric 
change of the plastic zone as a result of external load. This can be done simply by imagi- 
narily slicing the sample longitudinally along the loading direction and monitoring the 
two-dimensional change of the plastic zone in each plane. 

3. Experimental Results 

The discussion of the experimental results will be divided into four parts: (a) 
dislocation density vs. distance, (b) slip line analysis, (c) FEM results, and (d) voids. 

3.1. Dislocation density vs. distance 

3.1.1. 20 V% SiClAl 
composite 

The plot of disloca¬ 
tion densities (p) 
versus distance from 
the fracture surface 
(x) is shown in Fig. 8. 
The value of p is 
higher near the 
fracture surface region 
and decreases as x (ie. 

Fig. 8 Dislocation density vs. distance from fracture surface for 
20% SiC/Al composite. 
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the distance from the fracture surface) increases. When the distance is 3 mm and beyond, p 
is constant and the value of p is approximately equal to the value of p of an undeformed 
sample. This value is about one order of magnitude lower than that at the 0.3 ~ 0.5 mm 
region, i.e. near the fracture surface region. 

In general, the values of p are higher in the whiskered and heat treated samples. 

The difference between the dislocation densities of each micrograph in the same foil 
increases as x increases, which means the dislocation densities are more uniformly 
distributeed near the fracture surface region. 

In some foils, the difference in the dislocation densities is quite large (may exceed two 
orders of magnitude). To avoid confusion, the reported data are the average values of 
micrographs taken from each location rather than those of a single micrograph, i.e. there 
will not be two orders of magnitude difference deviation shown in the plot. 

Several attempts were made to deform 20 V% SiC/Al composites until a "neck" 
occurred, then stopping the test prior to fracturing the sample, and then cutting TEM foils 
from the sample in the "neck" region. All attempts were failures; either the sample 

fractured before the 
test was stopped or 
upon examination of 
the TEM foils. There 
was no difference in 
density between the 
perceived "neck" 
region and away from 
the "neck". 

3.1.2. 5 V% SiC/AI 
composites The plot 

Fig. 9 Dislocation density vs. distance from fracture surface of dislocation 

for 5 V% SiC/Al composite. 














Fig. 10(b) A transmission electron microscopy micrograph of foil taken from the neck 
section of 5 V% whisker SiC/Al tensile sample. 
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Fig. 11 A plot of dislocation density in the neck region and 
the end sections of a tensile sample which has been deformed 
to produce a neck. 


densities (p) vs. 
distance from the 
fracture surface (x) is 
shown in Fig. 9. 
Similar to 20 V% 
(Fig. 4) at 0.3 ~0.5 
mm, i.e. near the 
fracture surface 
region, p is higher and 
approximately equals 
the 20 V% samples at 
the same region, then 
p decreases as x 
increases. But 
contrary to 20 V%, at 


3 mm and beyond, p rather than staying constant and the values of p are 3 to -5 times 
higher than that of undeformed samples. 

Also, in agreement with 20 V%, p is higher in whiskered and heat treated samples. 

As in the case of 20 V% SiC/Al composite, several samples were tested until a "neck" 
had formed. Then, TEM foils were cut from the "neck" region and the ends of the samples. 
Figure 10(a) is a TEM micrograph of the end region of the sample and Fig. 10(b) is a TEM 
micrograph showing the high dislocation density in the "neck" region. The variation of the 
dislocation density from the "neck" region to the ends of the sample is shown in Fig. 11. 

3.1.3. 0 V% SiC/Al samples 

The plot of dislocation densities p vs. distance from the fracture surface (x) is shown 
in Fig. 12. Near the fracture surface region (0.3 to - 0.5 mm), p is higher, but the value of 
p is about one third of 20 V% samples at the same region and unlike the 20 and 5 V% 
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samples, as x 
increases p decrease.* 
in a much slower rate. 
The values of p in this 
plot are about 2 
orders of magnitude 
higher than an 
undeformed sample. 

Similar to 20 and 5 
V% samples, p is 
higher in heat treated 
samples. 

Fig. 12 Dislocation density vs. distance from fracture surface 3.2 Slip line analysis 

for 0 V% SiC/Al composite. 

The development of 
slip lines, i.e. the 

initial slip lines, began in a region which contained a higher than average density of SiC 
particles (that is, a cluster). The very first slip lines which occur in this higher density of SiC 
particle region, occur at the corners of SiC particles, as shown in Fig. 13. Figures 14,15, 16 
are micrographs taken at successive loadings. In other areas of the sample there are no or 
very few slip lines. 

3.3. FEM 

3.3.1. Two-Dimensional FEM 

The global stress-strain curves of the composites produced by calculation are shown in 
Fig. 17. It is clear that the more the whiskers are clustered together, i.e. the smaller the Sy 
and Sh in Fig. 6b, the lower the yield stress. This is consistent with the experimental data 
described from previous section; that is, plastic deformation is initiated within the region 
where there is a higher localized particle concentration (cluster). If a periodic arrangement 
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STRESS-STRAIN RELATION 



Fig. 17 Stress strain 
curves produced as a 
result of clustering 
where curve 1: Sv/2V 
= 1, V 2H = 1 ; 
Curve 2: S v /2V = 
0.75, S h /2H = 0.429, 
Curve 3: Sy/2V = 0.5. 
S h /2H = 0.143; 

Curve 3: Sy/2V = 0.5, 
S h /2H = 0.143; Curve 
4: S v /2V = 0.25, 
S h /2H = 0.071. 


of the reinforcement particles represents a more uniform plastic deformation best 
achievable by the composites, clustering of these particles would lead to nonuniformity of 
plastic deformation, that is, highly localized dislocation generation and, in this case, motion 
within the cluster of whiskers. It would, therefore, result in reduction of composite flow 
stress. 

3.3.2. Three-Dimensional FEM 

It was found that the plastic zone extended about one particle diameter from the edge 
of the particle which is in agreement with the experimental results [23]. Also, it was found 
that the plastic zone expanded transversely upon loading the sample in tension which is in 
agreement with the slip line analysis, whereas longitudinally it contracts near the whisker 
ends due to the compressive nature [24] of the plastic zone at that part of the matrix. 
However, this process generally means the generation of dislocations of opposite sign; 
annihilation of dislocations does not make a significant contribution to the process. 
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3.4. Voids 

The TEM micrographs were examined to determine the density of 
voids at the interface between the matrix and the SiC whiskers or particle. Voids were 
found, however, the density of voids near the fracture surface was greater than the density 
of voids in the non-deformed regions of the samples. 

4 . Discussion and Conclusions 

The discussion of the experimental results will be related almost exclusively to 
explaining how and why the deformation is so localized in the SiC/Al composites. There 
will be a brief discussion of the lack of voids. 

The evidence obtained from the dislocation density and the slip line investigation is 
overwhelming: deformation in higher volume fraction composites is very localized. From 
an analysis of the configuration of SiC particles in the samples which were used for the slip 
line analysis, it is evident that slip lines began in regions which contained a clustering of SiC 
particles. The slip propagated from one cluster region to another and deformation in the 
remainder of the samples is almost zero. The phenomena of localized deformation was 
duplicated in a FEM analysis. 

Although, by considering the fact that plastic deformation initiates within the cluster, 
the drop in yield stress can be easily explained. The question still remains to be addressed: 
Why should plastic deformation initiate in the clustered region where there is a local high 
volume fraction of reinforcement whereas it is obvious the higher the volume fraction of 
reinforcement, the stronger the material and the less prone to plastic deformation? 

To fully explain this phenomenon, it is necessary to consider the fact that deformation 
compatibility has to be satisfied when there is an inhomogeneous distribution of reinforce¬ 
ment particles, i.e. clustering. Since the clusters are stiffer, upon satisfying compatibility, 
the stress would be distributed in such a way that the clusters would bear more load than the 
rest of the matrix. Secondly, since the local volume fraction is higher within clusters, for the 
same amount of deformation, it must activate more slip systems in the matrix to accoinmo- 
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Fig. 18(a) This shows localized effective plastic strain at 
various locations vs. global strain. "Homogeneous material 
applies to a homogeneous material subject to a tensile load. 


date the same amount 
of deformation. 
Figure 18(a) shows 
effective plastic strain 
at various locations 
within the matrix [Fig. 
18(b)] as a function of 
global strain. If a 
sample of homoge¬ 
neous material was 
tested in tension, then 
the local strain would 
be the same as the 
global strain. This 
behavior is indicated 
by the line labeled 


Fig. 18(b) Locations 
where effective plastic 
strain is taken. 
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Fig. 18)c) Effective plastic strain vs. global strain at whisker 
ends for clustered and regular arrays (at location 1) for 
clustered array as shown in Fig. 18(b). It is clearly indicated 
that the clustered array has a higher effective plastic strain rate 
for the same global strain increment, which simply means 
clustered arrangement of the reinforcement particles is more 
susceptible to external loading in terms of plastic deformation. 


homogeneous 
material" in Fig. 18(a). 
From the plot, it can 
be easily seen that 
plastic strain rate 
between the fiber 
ends in the cluster is 
much higher than 
elsewhere in the 
matrix. Fig. 18(c), 
shows a comparison of 
the regular and a 
cluster whisker dis¬ 
tribution. The regu¬ 
lar distribution also 
has a higher plastic 
deformation rate at 
the whisker ends, but 
the cluster distribu¬ 


tion has a higher rate. This signifies that each cluster is being treated as a reinforcement 
particle since apparently it has a higher it has a higher Young’s modulus, which is in 
agreement with earlier work [24], and the stress level in the matrix within the cluster. 
Therefore, more plastic deformation will be induced. 

The lack of voids at the matrix/whisker interface is probably related to several 
phenomena. First, it has been shown [25] that the bond between SiC and A1 is very good. 
This bond is several times stronger than the bond between Fe 3 C and Fe. Second, there is a 
region of cold work, i.e. a high dislocation density region, around the whisker. As a result of 
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this hardened region around the whisker, a triaxial stress develops in the matrix some 
distance away from the whisker. It is into this region that the crack propagates. The end 
result is that very few voids are found at the interface between SiC and Al. 
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STRENGTHENING OF COMPOSITES DUE TO MICROSTRUCTURAL CHANGES 

IN THE MATRIX* 

By 

R.J. Arsenault, L. Wang and C.R. Feng 
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ABSTRACT 

The addition of discontinuous silicon carbide (SiC) to aluminum (Al) alloys can result in 
a five-fold increase in the yield stress. The magnitude of the increase is obviously a function 
of the volume fraction and the particle size of the SiC. Previously, it was proposed that the 
strength increase due to SiC addition to Al alloys was the result of change in the matrix 
strength, ie. an increase in dislocation density and a reduction of subgrain size. The data 
obtained from a series of experiments indicate that dislocation density increases with an 
increase in volume fraction of SiC and decreases with an increase in particle size. The 
subgrain size decreases as the volume fraction increases and increases as the particle size 
increases. 

There is a good correlation between the microstructural changes in the matrix and the 
changes in the yield stress of the composites. 

1. INTRODUCTION 

The concept of strengthening due to an increase in dislocation density has been in 
existence for several decades, probably since shortly after the original proposal of edge 
dislocations by Orwan, Polyani and Taylor [1]. Over the subsequent years, there have been 
numerous investigations which have correlated a strength increase with an increase in 
dislocation density. The strengthening due to a reduction in subgrain size is not as old as 
the concept of strengthening due to an increase in the dislocation density [2]. However, the 
correlation has been demonstrated quite clearly by several individuals; in particular, by 


* This research was supported by the Office of Naval Research under contract No. 
N00014-85-K-0007. 












McQueen et al. [3]. 

It had been proposed by Arsenault [4] a few years ago that the strengthening of Al and 
Al alloys due to the addition of discontinuous SiC reinforcements occurred as a result of 
generation of dislocations which produced an increase in dislocation density and a reduced 
subgrain size in the matrix. The generation is due to a difference in thermal coefficient of 
expansion between the matrix and the reinforcement (aCTE). The presence of the high 
dislocation density within the matrix has been demonstrated along with a very small sub¬ 
grain size. An in situ investigation has demonstrated the validity of the aCTE mechanisms 
as a means by which the high dislocation density is produced within the matrix [5]. 

Further, it has been shown that a very simple model based on prismatic punching is 
capable of producing an adequate dislocation density [6]. However, it should be clearly 
pointed out that this model will only predict the lower limit of the dislocation density. 
Attempts to calculate the upper limit of the dislocation density due to the aCTE effect have 
proven unsuccessful, for all of the calculations predicted a very high dislocation density, i.e. 
approaching infinity, which obviously is not correct. 

However, there remain several questions pertaining to the details of the strengthening 
mechanism of discontinuously reinforced metal matrix composites based on the strengthen¬ 
ing of the matrix due to a change in the microstructure of the matrix. 

The first question is related to transmission electron microscopy (TEM) procedures 
for determining dislocation densities; a) Does the sample preparation technique, i.e. 
dimpling and ion milling, introduce damage into the TEM foil? and b) If a thin region of 
the foil is examined (i.e. the viewing region) at 100-200 KV, will a lower dislocation density 
be observed? There have been a few investigations which have determined the dislocation 
density as a function of foil thickness and the most recent are by Fajita et al. [7-9]. In the 
case of pure Al, it was found that the foil thickness had to be greater than 1 urn to obtain a 
bulk dislocation density. The question arises, since SiC/Al is not high purity Al, does this 
requirement of foil thickness still apply? 









The second question: Is there an overall relationship between dislocation density (p), 
subgrain size ), reinforcement particle size (D) and volume fraction (V), and the increase 
in the yield stress of the composite? There have been several investigations where the 
dislocation density was determined for a specific particle size and volume fraction, but there 
have been no systemic investigations of these parameters. Also, there are no systematic 
investigations of the subgrain size as a function of particle size and volume fraction. 

Another question: Is it possible to increase the strength of the matrix sufficiently to 
account for the strengthening? In other words, if the A1 matrix alone (with no SiC present) 
had the same dislocation density as in the composite, would the A1 matrix have the same 
strength as the composite? Is it possible to increase the dislocation density and reduce the 
subgrain size within the A1 alloy matrix by cold working to approach the yield strength of a 
composite containing 20 V% reinforcement in the same A1 alloy matrix. 

An investigation was undertaken to address the questions raised above. 

2. MATERIALS AND TESTING 

A discussion of materials and testing will be divided into three portions corresponding 
to: 

2.1 Composite Materials 

2.2 Cold Rolled 0 V% 1100 A1 Alloy 

2.3 High Parity A1 
2.1. Composite Materials 

All of the composite materials were produced by a powder metallurgical procedure. 
The A1 alloy powders were mixed with particulate or whisker SiC, hot-compacted and 
hot-extruded into 12.7 mm diameter rods or hot pressed into a plate. The 0, 5 and 20 V% 
SiC whiskers/6061 A1 alloy composites plus the 0 V%/1100 A1 alloy composites were 
purchased from ARCO Silag*. The 20 V% SiC-250 pm particulate/1100 A1 alloy composite 


* Now, Advanced Composite Materials of Greer, South Carolina. 










was purchased from DWA* in the form of a plate. The remainder of the material was 
produced in the Metallurgical Materials Laboratory of the University of Maryland in the 
form of an extruded rod. 

From the extruded rod and plate, tensile samples were machined into a configuration 
as described elsewhere [4]. Prior to tensile testing, the samples were annealed for 12 hours 
at 803 K and furnace cooled. All of the tensile testing was conducted at room temperature 
at a strain rate of 10* 4 sec* 1 . 

Short sections of extruded rod (12 mm) and portions of as-pressed plate were 
annealed at 803 K for 12 hours and furnace cooled. From these sections and portions, 0.5 
mm thick disks were electrically discharge machined. The procedure for thinning the TEM 
foils and the method of determining the dislocation densities is given elsewhere [5]. The 
electron microscopy investigation involved the use of two microscopes, a JEOL CFX1 100 
operating at 100 KV and the Argonne National Laboratory 1 MEV microscope operated at 
1 MV. 

2.2 Cold Rolled 0 V% 1100 A1 Alloy 

The cold rolled strips were produced by taking 0 V% 1100 A1 alloy material pur¬ 
chased from ARCO Silag and rolling without any intermediate annealing until a reduction 
in thickness of ~ 90% was obtained. From the cold rolled strips, "dog bone" shaped tensile 
samples were machined with a 24.5 mm gage length and a width of 2.6 mm and a thickness 
of 1.42 mm. 

The tensile sample preparation and testing were the same as described in Part A. 
This also applies to the TEM foil preparation and examination. 

2.3. High Purity A1 

Relatively high purity A1 (99.99%) was obtained in the form of a rod, 12.7 mm dia. 
from ALCOA as a gift. This material was annealed for 12 hours at 803 K and furnace 
cooled. The TEM foils of high purity A1 were produced by two different procedures. Some 

* Composite Specialties of Chatsworth, California. 










133 




of the foils were produced using the electrical discharge machine, dimpling and ion milling 
procedure, using the same operating parameters as used for the composite TEM foils 
(Part 2.1). The other samples, after cutting into 0.5 mm thick disks by electrical discharge 
machining, were electrochemically thinned using a jet polishing device. The foils were 
examined at 100 KV and 1 MV. 

3. Results 

The discussion of the results will be divided into four portions: 

3.1 Yield Strength vs. Particle Size and Yield Strength vs. Volume Fraction 

3.2 Dislocation Density Determination 

3.2.1. Dislocation Density Due to Foil Preparation Procedures and TEM Operating 
Voltages 

3.2.2. Dislocation Density and Subgrain Size vs. Particle Size and Volume Fracture 

3.3. Changes in Dislocation Density and Subgrain Size due to Cold Rolling and the 
Determination of the Yield Strength 

3.1 Yield Strength 

Of interest is the incremental increase in the yield stress (which is defined as the stress 
at a plastic strain of 0.2%) due to the addition of the reinforcement. Therefore, the 
yield stress of the 0 V% matrix alloy is subtracted from the yield stress of the composite. 
The incremental yield stress (Aoy) is what is plotted and listed in all cases. The variations 
in Aoy of the SiC whisker/6061 A1 alloy composites were considerably less than in the 
SiC/1100 A1 alloy composites. This is simply due to the greater homogeneity of the rein¬ 
forcement within the matrix in the 6061 A1 alloy composites. 

As expected, as the size of the SiC particulate increased, the strength decreased, as 
shown in Fig. 1. The squares represent the average of at least three tests and there could be 
as many as ten tests for a given particle size. The arrow bars define the range of scatter. 
The ultimate strength is only slightly larger than the yield strength indicating that the work 
hardening rate past yielding is small. The uniform strain was fairly constant, between 3 and 
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4% for the 0.5 to 9 composite materials, then increasing to 6-7% uniform strength for 
the 20 and 70 fim particulate and slightly larger on average for the 250 /im particulate size 
composites. 

Increasing the volume fraction of SiC resulted in an increase in strength (Fig. 2), but 
not linearly as predicted by the classical-continuum-load transfer-shear lag model [10], 
which is indicated by the dashed line in Fig. 2. 

3.2 Dislocation Density Determinations 

3.2.1. Dislocation Density Due to Foil Preparation Procedures and TEM Operating 
Voltages 

The two different series of investigations were undertaken to determine if the TEM 
foil preparation procedure was introducing damage into the composite. The first investiga¬ 
tion involved preparing TEM foils from annealed and furnace cooled high purity A1 rods by 
two difference routes. The first involved electrical discharge machining, followed by 
dimpling using 0.25 qm diamond paste and then followed by ion milling, i.e. the same 
procedure which is used in the composite. The second technique involved electrical 
discharge machining and then followed by electro-chemically thinning. Figure 3(a) is a 
micrograph taken of an ion milled thinned sample and Fig. 3(b) is taken from a sample 
which was electrochemically thinned. As can be seen, the dislocation densities are low in 
both cases. The only recognizable difference between the ion milled and the e’ectro- 
chemically thinning is the "mottled" appearance of the ion milled micrograph. This mottled 
appearance is the result of the non-uniform thickness of the ion thinned foil; the reason for 
the non-uniformity is not known. 

However, it is safe to conclude that the dimpling and ion milling do not introduce any 
detectable increase in the dislocation density. 

The second investigation was related to foil thickness (viewing thickness) as related to 
dislocation density. If the operating voltage of the TEM is 100 KV, the thickness of the foil 
in which dislocations may be viewed is 0.2 to 0.4 ^m. However, if the operating voltage is 










1 MV, the thickness of the foil increases to 1 to 2 pm. If we consider the micrograph shown 
in Fig. 4, we can see the effect of foil thickness on the dislocation densities. Figures 4(a) 
and (b) are micrographs taken at 100 KV and the foil thickness in Fig. 4(b) in which a few 
dislocations are visible is 0.3 pm and the dislocation density is 1.6X10 13 m' 2 . Fig. 4(c) and 
(d) are micrographs taken at an operating voltage of 1 MV. The dislocation density in Fig. 
4(d) is 4.5X10 13 m' 2 and the foil thickness is 0.9 pm. Figures 4 (a), (b), (c) and (d) are 
micrographs taken from the same area of a foil under the same operating diffraction 
condition, only a thinner area around the brim of the hole is transparent for the 100 KV 
microscope. This is a general result for all composites examined. If a low operating voltage 
is used, i.e. 100-200 KV, the apparent dislocation density is low compared to the dislocation 
density obtained from foils examined at 1 MV. 

Therefore, in order to obtain realistic measures of dislocation density, it is necessary 
to examine thickness foils which needs a operating voltage (800 KV to 1 MV). 

3.2.2. Dislocation Density and Subgrain Size vs. Particle Size and Volume Fracture 

The determination of dislocation density by TEM is fraught with many perils. Two of 
the many possibilities have already been addressed, i.e. introduction of dislocations during 
sample preparation and the effective viewing of the foil thickness on the observed density. 
There are certainly several others, especially if numerous slip systems are operative as they 
certainly are in composites. In a two-beam condition, this means that several operative 
g-vectors have to be considered, but when operating at 1 MV, a many-beamed condition can 
easily be used and a larger fraction of the dislocations present become visible. However, it 
is still necessary to tilt the sample and take several different pictures. Also, there is an 
inherent variation of dislocations from place to place within the matrix. For example, the 
dislocation density is higher near the particles and decreases with distance away from the 
particle. This variation is especially critical in large particle size and low volume fraction 
samples. In order to overcome this difficulty, it is necessary to take measurements of many 
micrographs. The procedure for measuring the dislocation density is given elsewhere [5]. 
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The data given for each condition (i.e. particle size or volume fraction) in Figs. 5 and 6 
represent tens to hundreds of micrographs of a given foil and 5-6 different foils were 
examined for each condition. The dislocation densities reported are the average dislocation 
densities in the matrix of annealed samples. As stated above, the dislocation density is not 
uniform in the matrix. 

Upon increasing the volume fraction, the dislocation density increases, but upon 
increasing the volume fraction, the subgrain size decreases as shown in Fig. 5, and upon 
increasing the particle size, the opposite is the case, the dislocation density decreases (and 
the subgrain size increases) as the particle size increases (Fig. 6). 

3 3 Effect of Cold Rolling 

The yield strength and the dislocation density were determined of a 0 V% 1100 A1 
alloy which had been cold rolled ~ 90%. The cold rolling reduced the subgrain size 
slightly from 5-10 nm to 2 ^m, but there was a significant increase in dislocation from less 
than 10 1 °m to 1.5X10 13 m. The cold rolling increased the yield strength from 30 MPa to 
146 MPa. If a comparison is made with the 9 nm SiC particle size to 20 V% 1100 A1 alloy 
matrix composite, it can be seen that the Aoy is greater for the cold rolling material than 
that of the SiCp/1100 AI alloy composite (Table I). However, the interesting point is that 
the dislocation density in the composite is much greater than that of the 0 V% 1100 Al alloy 
cold rolled material. 


4. DISCUSSION 


The discussion will follow along the lines of the questions raised in the Introduction. 
4.1 Dislocation and Subgrain Size Determinations 

From the investigation of high purity Al, it was shown that the dislocation density in 
the foils, produced by dimpling and ion milling was the same as that in the electro- 
chemically produced foils. A main concern of ion milling is that the heat generation occurs 
and that this heat generation would result in a decrease in the dislocation density. In terms 
of dimpling, which involves grinding the foil with 0.25 nm diamond paste, this could 
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plastically deform the sample, i.e. introduce dislocations. These dislocations would have an 
entirely different appearance as compared to the remnant dislocations which remain from 
the long term high temperature annealing. The outcome of this investigation is thaf^the 
dimpling and ion milling do NOT introduce any dislocations into the TEM foil. 

If the foil thickness increases, the density of dislocations should increase up to 
center foil thickness and remain constant. For if we consider a homogeneous material, for 
example, high purity AI. The reason for the decrease in the density in a thinner foil is that 
the presence of the image force should move a proportionally larger number of dislocations 
from the thinner foil. 

In the case of SiC/Al alloy composites, we are not dealing with a homogeneous 
material. There are thermal residual stresses within the sample prior to thinning. If we 
consider points labeled 2 in Fig. 7, there is a tensile residual stress acting in the direction 
C 1 , A 1 and B 1 . If the sample is cut along the lines A-A 1 and B-B 1 , the stress at the 
surfaces will have to relax to zero. A simple manner of stress relaxation is by dislocation 
motion, i.e. dislocation motion out of the foil. Again, the thinner the foil, the more exten¬ 
sive the proportional dislocation loss out of the sample. The image force is still present and 
would tend to reduce the dislocation density still in the thinner area of the foil, just as in the 
homogeneous material. Therefore, in order to obtain reasonable measures of the disloca¬ 
tion density, it is necessary to examine thicker portions of the TEM foil, which means 
examining the foil at 1 MEV. The determination of the subgrain size is not as dependent on 
foil thickness or the foil preparation procedure. 

4.2 CORRELATION BETWEEN MICROSTRUCTURAL CHANGES AND YIELD 

STRESS 

The increase in the yield stress with an increase in the volume fraction or the decrease 
of yield stress with an increase in the particle size is to be expected. Correspondingly, the 
dislocation density and the subgrain change in the expected manner. However, it is neces¬ 
sary to consider if THE CHANGES IN THE YIELD STRESS CAN BE RELATED TO THE 








CHANGES IN DISLOCATION DENSITY AND THE SUBGRAIN SIZE. 


Before considering these correlations, it is necessary to point out again that disloca¬ 
tion density within the matrix is not uniform, and secondly, there can be a significant 
difference in stress at the proportional limit and the yield stress defined at a plastic strain of 
0.2%. This would indicate that there is significant work hardening on a local scale occur¬ 
ring within the matrix. The reason for stating that it is on a local scale is due to the fact 
that very little plastic strain is occurring between the proportional limit and the defined 
yield stress. Investigations have shown that it is very difficult to determine any changes in 
the overall dislocation density at yielding [11]. It is not known at the present time where the 
predominance of dislocation motion has occurred when a tensile plastic strain of 0.2% has 
been obtained, a FEM analysis [11] seems to indicate that the predominance occurs in the 
higher dislocation density region. Therefore, a assumption is made that the averag e 
increase in the dislocation density is small due to the initial loading to a plastic strain of 
0.2%, and that the stress field which opposes dislocation motion had a strain of 0.2% which 
is due to the average dislocation density of the annealed undeformed sample. Therefore, 
we can use the well-established relationship 

AOy = 

where Ao y is the increase in the yield stress of the composite over that of the 0 V% matrix 
material, « is a constant and Hansen [12] has shown that for Al, this 1.25, p is the shear 
modulus of the matrix (2.64x 10-4 MPa), b is the Burgers vector (2.86x 10" 10 m) and p is the 
dislocation density. 

Then we can calculate the incremental increase in the yield stress due to the presence 
of the increase in dislocation density. This is the increase which should be added to the 
annealed 0 V% Al alloy matrix yield stress. In Table II are listed the incremental increases 
in the strength due to an increase in dislocation density as a function of changes in particle 
size at a constant volume fraction of 20%. 
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If we consider the increase in dislocation density with the increase in volume fraction, 
it is obvious that the strengthening due to the increase in dislocation density alone is very 
large (Table III). In this case, it is necessary to subtract out the dislocation density of the 
0 V% material, for this is appreciable when considering smaller volume fractions. In a case 
of particle size increase, this subtraction is not necessary, for the dislocation density of the 
0V% annealed 1100 A1 alloy matrix was very low. 

In terms of strengthening due to a reduction in subgrain size, it is necessary to rely on 
the data of McQueen. In Table II are listed the predicted increase in strength due to a 
reduction in the subgrain size AogQ. The subgrain size of the 1100 A1 matrix is much larger 
than the subgrain size of the composite so a subtraction is not necessary. The amount of 
strengthening due to a reduction in subgrain size in the small particle size range is greater 
than that in the large particle size range. If we consider the effect of volume fraction on the 
change in the subgrain size, in this case, it is necessary to subtract the increment of 0 V%. 
The increase in strength due to the reduction in subgrain size becomes predominant at 
20 V%. 

In general, the predicted increase in the yield strength is much larger than the 
experimentally determined yield strength. There are several possible explanations for this 
discrepancy; a) the value taken for the constant o in Eq. 1 is too large. The values 
reported in the literature for a range from 0.5 to 1.25. The second possibility is that the 
assumption of using an average dislocation density is not valid for the plastic strain at 0.2% 
and may be due to the motion of dislocations in the matrix where the dislocation density is 
much lower than the average density. This certainly could account for the discrepancy in 
the 250 particle size composite and the 1 V% composite where the variation in the 
dislocation density is the greatest. 


43 STRENGTHENING OF THE MATRIX 


As discussed in the Introduction, there is a question as to whether the matrix can be 
strengthened to the required levels by changes in the microstructure, i.e. can there be a 
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sufficient increase in the dislocation density or a sufficient decrease in the subgrain size? 

The cold rolled investigation has shown that the 0 V% 1100 A1 alloy material can have 
a yield stress which is comparable in strength to a 20 V% SiC/Al alloy composite. 

5. CONCLUSIONS 


This investigation had as its purpose to address several questions concerning the 
observed strength of discontinuous SiC/Al alloys, especially related to the mechanism 
proposed which is based on the changes in the microstructure of the matrix due to differ¬ 
ences in thermal coefficient of expansion between the reinforcement and the matrix. 

From a consideration of the data, several conclusions were obtained: 

• The observed strength of the discontinuous SiC/Al composite can be accounted for in 
terms of a mechanism based on a change in the microstructure of the matrix. 

• This change in the microstructure, i.e. an increased dislocation density and a reduced 
subgrain size as compared to 0 V% matrix alloy, is due to the difference in coefficient 
of expansion between the reinforcement (SiC) and the matrix (Al). 

• The procedure of TEM foil preparation of the composites which involves electrical 
discharge machining, dimpling and ion milling does not introduce dislocations into the 
matrix. 

• In order to obtain realistic values of the dislocation density, it is necessary to view 
thick foils, i.e. an operating voltage of 1 MV. 
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TABLE I 

Comparison of 9x10‘ 6 m Particle Size 20 V% SiC/Al Composite 
with Cold Rolled 0 V% 1100 A1 Composite 


Material 

Dislocation 

Subgrain Size 

A Experiment 

A Calculate 


Density 


Yield Stress 

Yield Stress 


pxlO 14 m 2 

>xl0‘ 6 m 

Ao yo. 2 % MPa 

MPa 


1.1 

2 

88 

168 

0 V% CR 

.15 

2 

116 

106 


a calculated yield stress = the predicted increase in strength due to an increase in p and a 
decrease in X. 

a experimental yield stress = the yield stress of the composite minus the yield stress of 0 
V% matrix. 


TABLE II 

Predicted Increase in Yield Strength of 20 V% SiC 1100 Matrix 
and Experimentally Measured Values 


Particle Size 

D(/im) 

Aa p 

MPa 

A°SG 

MPa 

Aa ypred 

MPa 

A o v 
>exp 

MPa 

0.5 

126.6 

131 

257.6 

153 

9 

99 

69 

168 

88 

70 

79 

20 

20 

69 

250 

59.7 

0-13.8 

59.7-73.5 

26 


• 

Keys 

The yield stress defined at 0.2% plastic strain minus the 0.2% yield stress 

of the 0 V% 1100 matrix 


Aa Yexp 


Ao p 

The predicted increase due to the increased dislocation density 

• 

A(7 SG 

The predicted increase due to the reduced subgrain size 


A<7 y j 

*pred 

ThesumofAOp + Ao sg 
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TABLE III 

Predicted Increase in Yield Stress of 6061/AI Alloy 
as a Function of Volume Fraction of SiC Whisker 


Volume Fraction 

V% 

A o p 

MPa 

Ao S g 

MPa 

Ao ypred 

MPa 

A °yexp 

MPa 

1 

37.7 

0 

37.7 

6.9 

5 

56.6 

7.0 

63.6 

52 

20 

70 

48.3 

118.3 

100 


Keys 

A °yexp 

The yield stress defined at 0.2% 


the 0 V% 1100 matrix 

A °p 

The predicted increase due to the 

A°SG 

The predicted increase due to the 

A °y P red 

The sum of Ao^ x Ao sg 


plastic strain minus the 0.2% yield stress of 

increased dislocation density 
reduced subgrain size 
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Fig. 1 The change in the yield stress as a function of SiC particulate size. The matrix is a 
1100 A1 alloy. Tne composite is in the annealed condition. 
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Fig. 2 The change in the yield stress as a function of the volume fraction of SiC (SiC 
particulate is ~ 5 vm in dia.). The matrix is a 6061 A1 alloy. The composite is in the 
annealed condition. 
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Fig. 3(a) A transmission electron micrograph of high purity Al in the annealed condition 
thinned by dimpling and ion milling. 

















Fig. 4(a) A transmission electron micrographs taken at an operating voltage of 100 KV. 
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Fig. 5 The change in dislocation density as a function of volume fraction of SiC (SiC 
whiskers have a length to diameter on average of ~ 2 and the dia.~ 0.5 /im). 
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Fig. 6 The change in dislocation density as a function of particulate size. 










